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1.  Abstract 


A  reliable,  efficient  and  secure  electric  power  system  is  necessary  for  the  operation  of  critical 
buildings  in  a  base  or  the  whole  base  itself.  This  is  also  applicable  for  deployed  force  in  forward 
bases,  which  have  to  be  put  into  service  quickly  and  reliably.  At  present,  there  is  a  need  to  design 
a  distributed  and  autonomous  subset  of  a  larger  grid  or  a  microgrid  to  increase  the  security  and 
reliability  of  electricity  supply.  The  objective  of  this  work  was  to  model  and  simulate  a 
specialized  microgrid  called  an  Intelligent  Distributed  Autonomous  Power  Systems  (ID APS), 
which  play  a  crucial  role  in  building  a  scalable  power  grid  that  facilitates  the  use  of  renewable 
energy  technologies.  Microgrid  device  models,  including  distributed  energy  sources  and  loads,  as 
well  as  their  control  algorithms,  were  developed.  Several  case  studies  were  simulated  to  evaluate 
the  operation  of  the  ID  APS  microgrid  during  parallel  and  islanded  operation  modes.  Simulation 
results  indicated  that  the  proposed  ID  APS  control  model  was  able  to:  (i)  perform  demand 
management  during  nonnal  operating  condition;  (ii)  island  the  microgrid  from  the  main  grid 
once  an  upstream  fault  is  detected;  (iii)  secure  critical  loads  and  shed  non-critical  loads  according 
to  the  given  priority  list  during  emergencies;  and  (iv)  resynchronize  the  microgrid  to  the  main 
grid  after  an  upstream  fault  is  cleared. 


2.  Objective 

In  response  to  the  Statement  of  Need  (SON)  NUMBER:  SISON-08-04:  Scalable  Power  Grids 
that  Facilitate  the  Use  of  Renewable  Energy  Technologies,  we  proposed  to  model  and  analyze 
the  operation  of  an  Intelligent  Distributed  Autonomous  Power  System  (ID APS),  which  provided 
the  opportunities  for  load  control  and  dispatch  of  distributed  energy  sources,  especially 
renewables.  This  effort  resulted  in  the  intelligent  distributed  autonomous  power  grid  that  could 
integrate  renewable  energy  technologies  and  minimize  reliance  on  external  energy  resources,  and 
thereby  reducing  fossil  fuel  consumption.  This  capability  would  facilitate  the  implementation  of 
renewable  energy  projects  and  enhance  energy  security  and  reliability  for  the  mission-critical 
parts  of  military  bases  and  campus-type  facilities. 

In  this  proof-of-concept,  the  building  blocks  of  an  ID  APS  microgrid  comprised  both  physical 
and  cyber  systems  as  illustrated  in  Fig.  1.  The  physical  system  consisted  of  hardware  elements, 
including  distribution  circuits,  power  electronic  devices,  DERs  —  which  included  both 
distributed  generators  and  storage  technologies  —  and  loads.  DERs  in  consideration  included 
solar  photovoltaics,  wind  turbines,  microturbines,  fuel  cells,  and  battery  storage  technologies. 

The  cyber  system  consisted  of  a  group  of  software  agents  that  acted  as  a  central  decision  support 
unit.  The  cyber  layer  communicated  with  its  associated  physical  layer  through  addressable 
communication  interfaces,  such  as  Internet  Protocol  (IP)  addresses. 
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During  normal  operating  conditions,  the  ID  APS  microgrid  runs  in  parallel  with  the  local  utility. 
At  the  same  time,  it  optimally  coordinates  internal  loads,  distributed  energy  resources  (DERs) 
including  generation  and  storage  devices  to  address  any  operational,  environmental,  economic, 
or  security  constraints.  In  the  event  of  an  upstream  outage,  the  microgrid’s  control  architecture 
was  designed  to  isolate  the  microgrid  from  the  local  utility.  In  such  a  situation,  the  microgrid 
performs  load  control  and  activates  its  internal  generators  to  secure  critical  loads  based  on  a 
given  prioritized  list. 

Specifically,  the  objective  of  this  research  was  to  build  a  simulation-based  model  of  an  ID  APS 
microgrid.  Specific  tasks  accomplished  included: 

Task  1:  Development  of  ID  APS  device  models; 

Task  2:  Development  of  local  control  algorithms  to  control  ID  APS  device  models; 

Task  3:  Development  of  an  ID  APS  energy  management  system; 

Task  4:  Development  of  addressable  IP-based  communication  interfaces  for  connecting  ID  APS 
device  models  and  ID  APS  energy  management  system; 

Task  5:  Simulation  and  evaluation  of  the  microgrid  in  both  parallel  and  islanded  operations. 

In  response  to  the  Statement  of  Need  (SON)  NUMBER:  SISON-08-04,  the  proposed  proof-of- 
concept  work  aimed  at  designing  and  evaluating  a  specialized  microgrid  to  enable  the  following 
key  characteristics: 

•  Intelligent:  In  an  ID  APS  microgrid,  all  DER  devices  and  loads  were  an  integral  part  of 
the  network  and  could  be  controlled  locally.  The  intelligent  features  of  an  ID  APS 
microgrid  implied  its  ability  to  perfonn  look-ahead  dispatch  and  demand  side 
management.  In  the  event  of  a  commercial  power  outage,  the  highest  priority  loads  were 
served,  while  non-critical  loads  were  shed.  This  depended  upon  the  given  prioritized  list, 
available  internal  generation  and  expected  outage  durations. 
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•  Distributed:  DERs  are  dispersed  in  nature.  The  ID  APS  microgrid  networked  and 
interconnected  various  DERs  through  the  integration  of  the  IP-based  control  model.  This 
thus  eliminated  single  point  failures  by  allowing  dynamic  response  of  distributed  control 
strategies,  thereby  increasing  robustness  of  the  local  electric  power  system. 

•  Autonomous:  The  ID  APS  microgrid  was  designed  to  disconnect  itself  from  the  local 
distribution  utility  and  operate  autonomously  once  an  outage  was  detected.  The  isolation 
of  the  ID  APS  microgrid  from  the  utility  prevented  any  possible  cascading  failure.  More 
importantly,  it  ensured  the  continuous  operation  of  mission-critical  facilities. 

•  Plug  &  Play  and  Scalable:  Each  physical  element  in  the  ID  APS  microgrid  was  assigned  a 
unique  IP  address.  This  allowed  all  DERs  and  loads  to  communicate  their  operational 
parameters  with  the  system  controller.  As  a  result,  the  current  configuration  of  the  ID  APS 
microgrid  could  always  be  updated  when  additional  generators  and/or  storage  devices 
were  added,  removed  or  relocated.  Furthermore,  since  each  ID  APS  microgrid  was 
modular  in  nature,  several  ID  APS  microgrids  could  be  interconnected  and  configured 
such  that  they  could  be  the  building  blocks  for  a  more  resilient  regional  electric  power 
system. 

The  criteria  for  success  included  the  demonstration  that  the  proposed  microgrid  could:  (i) 
perform  demand  management  during  normal  operating  condition;  (ii)  island  the  microgrid  from 
the  main  grid  once  an  upstream  fault  is  detected;  (iii)  secure  critical  loads  and  shed  non-critical 
loads  according  to  the  given  priority  list  during  emergencies;  and  (iv)  resynchronize  the 
microgrid  to  the  main  grid  after  an  upstream  fault  is  cleared. 

This  proof-of-concept  work  would  lay  a  solid  groundwork  to  gather  background  information  on 
developing  device  models  within  a  microgrid,  and  defining  its  distributed  control  environment.  It 
is  possible  to  extend  the  modeling  and  simulation  based  on  this  work  to  explore  the  feasibility  of 
microgrid  development  in  a  campus-type  facility.  The  developed  simulation  platform  with  some 
modifications  could  be  used  to  simulate  and  analyze  the  impact  of  deploying  certain  types  of 
distributed  generation  and  load  control  algorithms  within  a  microgrid.  Therefore,  this  makes  it 
possible  to  quantify  and  analyze  the  risks  and  the  impact  of  various  generation  and  load  options 
before  the  actual  field  implementation. 


3.  Background 


3.1  SERDP  Relevant 

Within  the  DoD  and  the  new  Anny  Energy  Strategy  for  Installations,  one  of  the  major  objectives 
is  to  decrease  dependence  on  fossil  fuels  and  increase  energy  security  http://army- 
energv.hqda.pentagon.mil/programs/plan.asp.  To  achieve  this  objective,  it  is  necessary  to 
diversify  DoD  current  use  of  the  local  electric  utility  by  integrating  various  types  of  distributed 
energy  sources,  including  renewable  energy  systems  such  as  wind,  solar,  and  other  advanced 
non-polluting  Distributed  Energy  Resource  (DER)  technologies  (e.g.,  fuel  cells  and 
microturbines). 
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At  present,  the  renewable  energy  integration  into  the  electric  power  grid  only  allows  parallel 
operation  with  no  capability  to  provide  power  to  critical  mission  facilities  during  a  commercial 
power  outage.  To  overcome  this  problem,  it  is  possible  to  network  these  power  systems  together 
in  a  scalable  intelligent  power  grid,  generally  known  as  a  microgrid.  Microgrids  have  emerged  as 
a  promising  new  means  to  integrate  various  DERs,  together  with  smart  load  control  strategies, 
with  the  aim  to  increase  the  resiliency,  security,  flexibility  and  efficiencies  for  the  mission- 
critical  parts  of  military  bases  and  certain  campus  facilities. 


3.2  Previous  Work 

Discussed  below  is  previous  and  related  work  related  to  microgrid  research,  development  of 
DER  models,  together  with  their  control  and  communication  architecture. 

3.2.1  Microgrid  R&D 

A  microgrid  comprises  the  interconnected  distributed  generation  (PV,  wind,  diesel, 
microturbines,  fuel  cells,  etc)  -  along  with  energy  storage  devices  (conventional  batteries, 
hydrogen  storage,  flywheels,  etc)  -  and  controllable  loads  at  low-voltage  distribution  levels. 

Such  systems  can  operate  in  parallel  with  the  local  utility  or  in  an  islanded  mode  during 
emergency  conditions.  Microgrid  is  one  of  the  key  technologies  recommended  by  policy  makers 
drafting  technology  roadmaps  for  electricity  delivery  in  many  countries,  including  the  United 
States  [1,  2],  the  European  Union  [3]  and  Japan  [4].  Subsequent  studies  have  all  promoted  the 
microgrid  concept  citing  the  increased  reliability  and  power  quality  it  provides  to  the  local 
utility. 

At  the  time  of  writing  this  report,  there  were  two  major  efforts  for  microgrid  development: 
CERTS/Sandia  Labs  Microgrid  Test  Bed  [5]  and  US  Anny  CERL/Sandia  Labs  Energy  Surety 
Project  [6].  The  former  developed  the  microgrid  test-bed  demonstration  with  American  Electric 
Power.  The  test-bed  comprised  three  low  voltage  feeders  at  480V  and  a  couple  of  distributed 
generators  [7].  The  work  was  very  comprehensive  and  covered  many  aspects  of  generator 
controls  and  parallel  operation  of  various  distributed  energy  sources.  The  latter  addressed  the 
Energy  Surety  Microgrid  project.  It  was  one  of  the  very  few  reported  work  that  focused  on 
microgrid  development  at  a  military  base  with  the  objective  of  creating  analytical  tools  and  a 
methodology  to  evaluate  the  impact  of  infrastructure  disruption  on  base  missions  [8]. 

The  main  challenge  identified  in  microgrid-related  research  was  how  to  intelligently  control  and 
coordinate  various  DERs  to  take  advantage  of  their  inherent  scalability  and  robustness  features. 
In  order  to  simulate,  analyze  and  evaluate  such  a  feature,  building  blocks  of  a  microgrid  were 
developed,  including  models  for  distributed  energy  sources  and  the  microgrid  control  algorithms 
based  on  multi-agent  technology. 

Previous  and  related  work  related  to  these  topics  is  discussed  below. 

3.2.2  Models  for  Distributed  Energy  Sources 

a)  Solar  PV  model 

Solar  Photovoltaic  (PV)  is  a  technology  that  converts  sunlight  directly  into  electrical  energy.  The 
output  is  direct  current.  The  major  components  of  a  PV  system  comprise  a  PV  array,  a  battery 
storage  unit,  an  inverter  and  a  charge  controller. 
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Previous  work  that  developed  solar  PV  models  included  the  following  publications  [4-12]. 
Modeling  techniques  of  PV  panels  were  classified  into  two  types:  numerical  techniques  and 
analytical  techniques.  The  former  required  iterative  processes  and  mathematical  tools  to  solve 
the  implicit  exponential  equation  associated  with  diode  and  photovoltaic  devices.  The  latter 
involved  more  simplified  and  approximated  processes;  yet  was  proven  to  be  accurate  without 
introducing  significant  errors.  A  very  good  numerical  technique  was  published  in  [9]  that 
predicted  performance  of  a  solar  photovoltaic  generator  in  various  operating  conditions.  Many 
subsequent  publications  had  also  followed  this  similar  method  [10,  11,  12,  13,  14].  Accurate 
analytical  methods  for  the  extraction  of  solar  cell  model  parameters  were  proposed  in  [15,  16]. 
Along  this  same  line,  Analytical  techniques  to  determine  parameters  of  solar  cells  were  proposed 
in  [17], 

b)  Wind  turbine  model 

A  wind  turbine  is  used  to  convert  power  in  the  wind  into  electricity.  Main  components  of  a  wind 
turbine  generator  are  turbine  blades  and  a  rotor,  a  gearbox  and  an  electrical  generator.  Turbine 
blades  convert  power  in  the  wind  into  mechanical  power  (torque)  that  is  transmitted  to  a 
generator  via  a  low-speed  shaft,  gearbox  and  a  high-speed  shaft.  The  output  of  a  wind  turbine 
can  be  a  DC  source,  an  AC  source  or  a  variable  frequency  AC  source,  depending  on  generator 
types.  A  power  converter  is  required  to  connect  a  wind  generator  to  a  utility  grid. 

The  modeling  and  control  of  wind  turbine  for  power  systems  were  described  in  many  previous 
publications  [13-16].  The  modeling  requirements  of  wind  turbine  for  a  power  system  were 
investigated  and  simulated  by  [18].  The  wind  turbine  model  and  its  perfonnance,  linearization 
and  control  were  presented  by  [19].  A  mathematical  and  dynamic  model  of  a  wind  turbine  with  a 
doubly-fed  induction  generator  for  grid  connection  was  presented  by  [20].  A  dynamic  model  of  a 
variable  speed  wind  turbine  with  permanent-magnet  synchronous  generator  to  be  connected  to 
the  grid  was  developed  in  Matlab/Simulink  by  [21].  By  integrating  the  state-of-the-art  modeling 
and  simulation,  Matlab/Simulink  offered  a  dynamic  wind  turbine  model  along  with  its  associated 
asynchronous  generator. 

c)  Microturbine  model 

Microturbines  are  small  combustion  turbines  that  produce  between  25kW  and  500kW  of  power. 
Microturbines  have  a  common  shaft  on  which  is  mounted  a  compressor,  a  turbine  and  a 
generator.  These  components  are  mounted  on  air  bearing,  thus  eliminating  friction  and 
maintenance  costs.  Microturbines  have  high  rotating  speeds  of  60,000  to  120,000  rpm. 

In  many  previous  works,  a  typical  microturbine  model  was  developed  by  connecting  a  gas 
turbine  model  to  a  pennanent  magnet  synchronous  generator  (PMSG)  model.  The  mathematical 
representations  of  a  heavy-duty  gas  turbine  model  were  proposed  by  [22].  The  model  was 
suitable  for  use  in  dynamic  power  system  studies  and  in  dynamic  analyses  of  connected 
equipment.  Based  on  Rowen’s  model,  authors  in  [23]  conducted  field  tests  to  derive  model 
parameters  and  validated  the  model’s  behaviors  with  actual  field  data  using  various  types  of  gas. 
Subsequently,  the  proposed  heavy-duty  combustion  turbine  model  by  Rowen  was  widely  used  as 
a  gas  turbine  section  of  the  microturbine  simulation  model  proposed  by  [24,  25,  26,  27].  In  these 
work,  the  circuit  models  of  the  permanent  magnet  synchronous  generator  were  developed  using 
the  rotor  reference  frame  (dq)  for  predicting  the  PMSG’s  transient  behavior. 

d)  Fuel  cell  model 
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Fuel  cells  are  electrochemical  energy  conversion  devices  that  convert  chemical  energy  in  fuel 
(H2  and  02)  into  electrical  energy.  By  products  of  fuel  cells  are  water  and  heat.  A  fuel  cell 
comprises  two  electrodes  (anode  and  cathode)  separated  by  an  electrolyte.  Fuel  (H2)  is  fed  to  the 
Anode.  Oxidant  (02)  is  fed  to  the  Cathode.  Power  is  produced  by  passing  of  ions  formed  at  one 
end  to  the  other  end  of  electrodes. 

Mathematical  and  dynamic  models  of  fuel  cells  were  described  in  many  papers  [46-49].  A 
dynamic  electrochemical  model  of  the  PEMFC  along  with  the  simulations  and  validation  results 
were  presented  by  [28]  that  showed  fuel  cell  dynamics  as  well  as  its  power  generation.  A 
mathematical  model  of  PEM  fuel  cells  was  presented  by  [29],  which  also  considered  the 
dynamic  response  of  the  fuel  cells  in  tenns  of  their  transient  and  thennal  properties.  A 
mathematical  model  of  PEM  fuel  cells,  together  with  their  control  architecture  and  a  power 
conditioning  system,  was  presented  by  [30],  A  dynamic  Simulink-based  model  with  validated 
results  in  both  steady  and  transient  states  was  presented  by  [31].  Furthermore,  Matlab  2008a 
already  implemented  generic  hydrogen  fuel  cell  stack  model,  which  represented  any  user-defined 
fuel  cells. 

e)  Battery  storage  model 

Batteries  are  used  to  store  excess  electrical  energy  from  the  generation  sources  and  supply  the 
load  during  the  times  when  the  sources  are  not  available.  Batteries  are  charged  when  an  external 
potential  is  applied  to  the  batteries’  terminals.  Batteries  are  discharged  when  an  external  load  is 
connected  to  the  batteries’  tenninals.  When  batteries  are  discharged,  the  batteries’  chemical 
reaction  is  reversed  and  the  absorbed  energy  is  delivered. 

Many  battery  models  were  developed  over  past  years  [32,  33].  A  Thevenin  equivalent 
mathematical  battery  model  was  developed  by  [34].  A  kinetic  battery  model  composed  of 
capacity  and  voltage  models  was  developed  by  [35]  which  represents  the  sensitivity  of  storage 
capacity  to  the  rate  of  discharge.  Several  other  battery  models  [36,  37]  were  developed  based  on 
the  physical  and  electrochemical  processes.  Based  on  some  of  these  previous  works,  a  complete 
battery  model  developed  by  [38]  available  in  Matlab/Simulink  comprised  a  controlled  voltage 
source  and  an  internal  resistance  with  current  discharge  characteristics. 

3.2.3  Multi-Agent  Systems  for  Microgrid  Control  and  Communications 

A  multi-agent  system  is  a  combination  of  several  agents  working  in  collaboration  in  pursuit  of 
accomplishing  their  assigned  tasks  to  achieve  the  overall  goal  of  the  system.  According  to  [39, 
40,  41],  a  software  program  was  declared  as  an  agent  if  it  exhibited  the  following  characteristics: 

•  Interacting  with  the  environment  and  other  agents 

•  Learning  from  its  environment 

•  Reacting  to  its  environment  in  a  timely  manner 

•  Taking  initiatives  to  achieve  its  goals,  and 

•  Accomplishing  tasks  on  behalf  of  its  user 

These  properties  signified  the  importance  of  a  multi-agent  system  in  developing  complex 
systems  that  enjoy  the  agent’s  properties  of  autonomy,  sociality,  reactivity  and  pro-activity  [42]. 
Multi-agent  systems  are  being  applied  in  many  systems  today.  Applications  of  agent-based 
systems  can  be  divided  into  single-agent  systems,  and  multi-agent  systems.  The  applications  of 
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the  former  include  situations  where  a  human  may  require  assistance  while  using  a  computer 
software.  These  are,  for  example,  infonnation  retrieval  and  filtering  software,  meeting  scheduler 
software,  mail  management  engine,  news  filtering  engine,  and  search  engine.  The  latter  is  where 
multiple  agents  work  together  to  achieve  a  particular  goal.  Examples  include  traffic  monitoring 
systems,  decision  support  systems,  manufacturing  systems,  telecommunications  and  network 
management  systems,  aircraft  maintenance,  military  logistics  planning,  and  power  systems. 

In  general,  multi-agent  systems  were  proposed  for  application  in  military  systems  [43,  44,  45], 
information  retrieval  systems  [46],  decision  support  systems  [47],  supply  chain  [48], 
transportation  [49],  communication  systems  [50]  and  many  more  [51,  52,  53,  54,  55].  In  the 
context  of  power  systems,  multi-agent  technologies  were  applied  in  a  variety  of  applications, 
such  as  to  perform  power  system  disturbance  diagnosis  [56],  fault  diagnosis  [57],  power  system 
restoration  [58],  power  system  secondary  voltage  control  [59]  and  power  system  visualization 
[60].  The  multi-agent  systems  for  power  engineering  applications  were  discussed  in  [61,62],  The 
distributed  control  approach  was  implemented  in  [63]  using  multi-agent  systems  technology.  In 
[64,  65,66],  the  authors  implemented  multi-agent  systems  for  optimal  operation  of  microgrids.  In 
[67],  a  multi-agent  system  was  proposed  that  attempts  to  restore  a  distribution  system  network 
after  a  fault. 


3.3  Proof  of  Concept  and  its  Technical  Challenges 

This  work  focused  on  modeling  and  simulation  of  a  microgrid  that  networked  variety  of 
distributed  energy  sources.  Key  technology  gaps  addressed  included  the  development  of  a 
distributed  agent  system  that  controlled  and  networked  distributed  energy  sources,  as  well  as 
analyzing  the  dynamic  response  of  distributed  control  strategies.  In  a  given  situation,  an  agent 
must  be  able  to  issue  a  control  signal  in  response  to  an  event  sensed  from  the  external 
environment  quickly  enough  to  manage  the  microgrid  in  a  timely  fashion.  The  simulation  and 
evaluation  case  studies  were  conducted  to  test  network  interoperability  in  both  on-  and  off-grid 
operation.  As  previously  discussed  in  Section  2,  our  success  criteria  were  to  demonstrate  that  the 
proposed  microgrid  could:  (i)  perform  demand  management  during  nonnal  operating  condition; 
(ii)  island  the  microgrid  from  the  main  grid  once  an  upstream  fault  is  detected;  (iii)  secure  critical 
loads  and  shed  non-critical  loads  according  to  the  given  priority  list  during  emergencies;  and  (iv) 
resynchronize  the  microgrid  to  the  main  grid  after  an  upstream  fault  is  cleared.  Through  the 
development  of  a  “plug-and-play”  interconnected  power  grid,  it  is  expected  that  this  work  can 
contribute  to  allowing  the  installation  or  deployed  force  to:  (1)  install  future  renewable  energy 
systems  and  (2)  effectively  control  and  optimally  benefit  from  the  power  that  is  generated  from 
DERs.  The  power  grid  should  provide  these  capabilities  during  both  normal  grid-connected  and 
emergency  or  islanded  operating  conditions. 
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4.  Materials  and  Methods 


4.1  Task  1:  Development  of  ID  APS  Device  Models 

A)  Task  Description 

Task  1  involved  the  development  of  ID  APS  hardware  models  using  the  SimPowerSystems 
toolbox  in  the  Matlab/Simulink  environment.  ID  APS  hardware  components  included  distribution 
circuit,  distributed  energy  resources  (DERs)  and  loads.  DERs  of  interest  included  solar 
photovoltaic  modules,  wind  turbines,  microturbines,  fuel  cells  and  battery  storage.  Below  are  the 
descriptions  of  how  the  model  of  each  physical  element  in  the  ID  APS  microgrid  was  developed. 

B)  Distribution  Circuit 

The  models  of  distribution  circuit  components,  including  circuit  breakers,  transfonners  and  lines, 
were  readily  available  in  the  SimPowerSystems  libraries.  These  existing  models  were  used, 
together  with  the  system  parameters  from  Virginia  Tech  Electric  Services  (VTES),  to  build  a 
simplified  distribution  circuit,  as  shown  in  Fig.  2.  System  parameters  needed  to  construct  the 
distribution  circuit  in  Matlab  included  system  voltage,  the  length  of  the  line,  line  resistance  and 
reactance,  transformer  size  and  the  circuit  topology. 


Fig.  2.  Distribution  circuit  of  interest 

C)  Distributed  Energy  Resources  (DERs) 


Description  of  the  development  of  each  model  is  summarized  below. 
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C.  1 )  Model  of  Solar  PV  Modules 

a)  Model  Description 

The  model  of  a  solar  photovoltaic  (PV)  generator  was  developed  in  the  Matlab/Simulink 
environment.  It  was  designed  such  that  standard  electrical  characteristics  of  solar  cells,  as  well  as 
temperature  and  solar  irradiation  were  used  as  inputs.  The  standard  electrical  characteristics  used 
as  inputs  to  the  model  included:  open-circuit  voltage  (Voc),  number  of  cells  connected  in  series 
(Ns),  short-circuit  current  (Isc),  number  of  cells  connected  in  parallel  (Np),  maximum  module 
power  (Pmax),  temperature  coefficient  of  open-circuit  voltage  and  temperature  coefficient  of 
short-circuit  current.  These  parameters  could  be  obtained  from  a  manufacturer’s  datasheet. 
Additional  inputs  to  the  model  were  cell  temperature  and  solar  irradiation,  which  could  be 
obtained  from  historical  data.  Model  outputs  were  PV  outputs,  including  voltage  and  current,  at 
the  maximum  power  point. 

b)  Mathematical  Formulae 

The  equivalent  circuit  of  a  photovoltaic  cell  is  represented  by  a  constant  current  source  (Iph) 
connected  with  a  diode  and  a  series  resistance  (Rs,cell),  as  shown  in  Fig.  3. 


Fig.  3.  Equivalent  circuit  of  a  solar  cell 

The  net  cell  current  (Icell)  can  be  determined  by  subtracting  the  diode  current  (Id)  from  the 
photo-generated  current  (Iph)  as  follows. 

i.,„  -  -h-  >„»  -  h  exP  ‘ ')  <Eci' 1  > 


where 

Icell  =  net  cell  current 
Vcell  =  cell  voltage 

Iph  =  photo-generated  current  which  depends  on  solar  irradiation 
10  =  diode  dark  satulation  current 
q  =  Magnitude  of  electron  charge  (1.6x10-19) 
k  =  Boltzmann’s  constant  (1.38x10-23) 

m  =  1-2,  a  constant  depending  on  material  and  physical  structure  of  a  solar  cell 
Tc  =  absolute  cell  temperature  (K) 

Let  Vt,Ceii  =  mkTc/q,  representing  a  cell  thermal  voltage.  With  the  assumption  that  the  photo¬ 
generated  current  (Iph)  equals  the  short-circuit  current  (Isc,ceii)  and  exp(q(V+IRs)/mkTc)  »  1 
under  all  working  conditions,  Eq.  1  can  be  rewritten  as: 
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(Eq.  2) 


I  cell  1 SC, cell  ^0  eXP 


^cell  +  hell^s,cell 


V, 


cell 


Considering  an  open-circuit  condition,  the  open-circuit  voltage  can  be  determined  from  Eq.  2 
given  the  solar  cell  generates  no  current  (Iceii  =  0)  as: 


^ OC,celI  ^ t, cell  M 


1  SC, cell 


(Eq.  3) 


As  a  result,  I0  can  be  expressed  as  a  function  of  Isc  as: 

(  V  ^ 

Wac^ex p-  oc’ce,,] 


\  ^ t  .cell  ] 

By  substituting  Io  from  Eq.  4  into  Eq.  2: 


I  cell  ^ SC , cell 


(  (v  -V  + 1  R 

1  v  cell  y  OC ,cell  x cellAVs, cell 

p  v 

\  \  v  t,cell  )) 


(Eq.  4) 


(Eq.  5) 


Eq.  5  represents  the  relationship  between  the  cell  current  and  the  cell  voltage  in  terms  of  short- 
circuit  current,  open-circuit  voltage,  as  well  as  series  cell  resistant  and  cell  temperature. 

A  PV  module  consists  of  many  solar  cells  connected  in  series  and  parallel  in  order  to  generate  a 
required  open-circuit  voltage  or  short-circuit  current.  Similarly,  a  PV  generator  consists  of  many 
PV  modules.  By  assuming  that  every  solar  cell  constituting  a  PV  module  or  a  PV  generator  are 
identical  and  has  the  same  operating  characteristic  under  the  same  illumination  and  temperature, 
the  characteristic  I-V  curve  for  a  PV  module  or  a  PV  generator  can  be  derived  by  considering 
that: 

I  Module  ~  I  cell  ’  ^  p  (^Q-  6) 

y Module  =  Kell  ’  N s  (Eq.  7) 


Where  iModuie  and  VModuie  respectively  are  the  current  and  voltage  of  a  PV  module.  Np  and  Ns  are 
the  number  of  cells  connected  in  parallel  and  series  that  make  up  a  PV  module.  IModuie  and 
VModuie  can  also  represent  the  current  and  voltage  of  a  PV  generator  when  Np  and  Ns  represent 
the  number  of  cells  connected  in  parallel  and  series  that  constitute  a  PV  generator. 


Using  Eq.  6  and  Eq.  7,  Eq.  5  can  be  rewritten  to  represent  the  relationship  between  the  module 
current  and  the  module  voltage  as: 


1  Module 


=  N  I 

iy  p1  SC, cell 


1  -  exp 


Module 


^ y OC, cell  +  I ModuteR-sjcellW  s  !  N p)^\ 


N.  ■  V. 


t,cell 


n 


(Eq.  8) 


This  is  the  key  equation  that  will  be  used  to  derive  the  I-V  curve  of  a  PV  generator.  Most 
parameters  on  the  right-hand  side  of  the  above  equation  ( V0c,ceii@STC ,  hc,ceii@STC  and  Pmax,ceii@STC ) 
are  readily  available  from  manufacturer’s  datasheet. 
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The  cell  series  resistance  (Rs,cell)  is  a  property  of  the  solar  cells,  and  must  be  determined  at  the 
standard  test  condition  (STC).  This  value  is  assumed  to  be  unaffected  by  temperatures  and  solar 
illumination.  The  cell  series  resistance  can  be  derived  using  the  following  relationships: 


R 


s,cell 


FF, 


STC 


FF  0 


STC  , 


X 


YoCji cell  @  STC 
I  SC, cell  @  STC 


(Eq.  9) 


where: 

FoC,cell@STC  Voc, module  /  Ns 
hc,cell@STC  Isc, module  /  Np 
P max,cell@STC  P max,  module  /  TVs'  Np) 

FFsTC  =  P max,cell@STC  /  (V OC,cell@STC  X  IsC,eell@STc) 

FFo,stc  =  (vocSTD  -  Ih(vocstd  +  0.  72))/(vocstd  +  1) 

Where  vocstd  is  normalized  voltage  at  the  standard  test  condition: 
VOCSTD  =  VoC,cell@STC  /  [k  (273+  Tc)  /  q] 


Derivation  ofVoc.ceii  and  Isc, ceil  at  an  Operating  Condition 

Some  adjustments  are  necessary  to  modify  open-circuit  voltage  and  short-circuit  current,  which 
are  given  at  the  standard  test  condition,  to  an  operating  condition  at  a  particular  irradiation  and 
temperature  level. 

It  is  generally  assumed  in  practice  that  the  short-circuit  current  of  a  solar  cell  depends  linearly 
and  exclusively  on  solar  irradiation.  In  addition,  the  short-circuit  current  of  a  solar  cell  increases 
with  the  cell’s  temperature.  Most  manufacturer’s  datasheet  gives  this  infonnation  in  terms  of  a 
short-circuit  temperature  coefficient  (A/°C). 


1  SC , cell  ^ SC  .cell  @  STC 


+  a(Tc-TsTC) 


(Eq.  10) 


Where: 

G  =  solar  irradiation  (W/m2) 

GSTC  =  solar  irradiation  at  the  standard  test  condition  (W/m2) 
a  =  temperature  coefficient  of  short-circuit  current  (A/°C) 

Tc  =  cell  temperature  (°C) 

TSTC  =  cell  temperature  at  the  standard  test  condition  (°C) 

The  open-circuit  voltage  of  a  solar  cell  decreases  with  the  cell’s  temperature.  Manufacturer’s 
datasheet  generally  gives  this  relationship  in  tenns  of  a  constant  C3  (V/°C).  In  addition,  the 
open-circuit  voltage  increases  slightly  with  irradiation  according  to  the  following  relationship. 
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inkT 


+  P(Tc-Tstc) 


(Eq.  11) 


^OC.cell  V OC,cell@STC  + 


-In 


q 


' STC  , 


Where: 

(3  =  Temperature  coefficient  of  open-circuit  voltage  (V/°C) 


C.2)  Model  of  Wind  Turbines 

a)  Model  Description: 

This  study  used  the  previously  developed  wind  turbine  model  in  Matlab  as  a  basis  for  the 
simulation  study.  In  this  case,  the  wind  turbine  generator  was  modeled  as  a  wind  turbine  section 
connected  with  an  electrical  generator  section.  The  turbine  section  was  modeled  based  on  the 
steady-state  power  characteristics  of  the  turbine  and  required  as  inputs:  (1)  the  generator  speed  in 
per  unit  of  the  nominal  speed  of  the  generator,  and  (2)  the  wind  speed  in  m/s.  Output  of  the 
turbine  block  was  the  mechanical  torque  (Tm)  in  per  unit  of  the  nominal  torque  of  the  generator 
that  applied  to  the  generator  shaft.  The  mechanical  torque  was  used  as  an  input  to  the 
asynchronous  generator.  The  output  of  the  wind  turbine  model  was  the  three-phase  power  that 
could  be  directly  connected  to  the  electrical  grid. 

b)  Mathematical  Formulae 

Eq.  12  describes  the  relationship  between  the  mechanical  power  output  (watts)  of  a  wind  turbine 
and  the  wind  speed  in  m/s. 

(Eq.  12) 


Where 

Pm  :  Mechanical  output  power  of  the  turbine  (Watts) 

Cp  :  Performance  coefficient  of  the  turbine 
|3  :  Blade  pitch  angle  (degree) 

A  :  Tip  speed  ratio  of  the  rotor  blade  tip  speed  to  wind 

p  :  Air  density  (kg/m3) 

A  :  Turbine  swept  area  (m2) 

v  :  Wind  speed  (m/s) 

Cp  represents  extractable  power  from  the  wind.  Cp  is  a  function  of  blade  pitch  angle  ( (  > )  and  tip 
speed  ratio  (X).  Maximum  possible  Cp  is  0.5928,  according  to  Betzz  limit.  In  this  study,  we  shall 
use  a  generic  equation  of  Cp  as  proposed  by  [68],  which  is  expressed  according  to  the  following 
formula. 


(A ,/3)  =  0.5176^ -0.4/3  -  5  j  •  e  A  +  0.0068A 


116 


•21 

aT 


(Eq.  13) 


Where: 


(  1  0.035  V 

{k  +  Om/3  13*  +  lj 
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According  to  the  above  equations,  the  maximum  possible  perfonnance  coefficient  of  a  turbine  is 
0.48  (Cpmax  =  0.48)  and  is  achieved  for  (3  =  0°  and  tip  speed  ratio  X  =  8.1.  This  value  will  be 
defined  as  the  nominal  value  (Xnom).  It  is  interesting  to  note  that,  at  a  constant  tip  speed  ratio, 
i.e.  X  =  8. 1,  increasing  blade  pitch  angle  will  result  in  lower  perfonnance  coefficient  (Cp),  thus 
reducing  the  turbine  power  output. 

Note  that  the  tip  speed  ratio  X  is  a  ratio  of  rotor  speed  to  wind  speed.  This  rotor  speed  is 
equivalent  to  the  generator  speed  because  the  turbine  rotor  and  the  generator  are  connected  via  a 
gearbox.  Eq.  13  implies  how  faster  the  blade  can  rotate  as  compare  to  the  wind  speed. 


2  nR^ 
A  = - 60_ 


v 


wind 


(Eq.  14) 


where 

R  :  Blade  radius  (meters),  i.e.  40m  for  a  2MW  V80  turbine 
vwind  :  Wind  speed  (m/s) 

rpm  :  Blade  speed  (rpm),  i.e.  9-19rpm  for  a  2MW  V80  turbine 

The  detail  of  asynchronous  generator  model  was  well  documented  in  Matlab/Simulink  and 
would  not  be  discussed  in  further  detail  here. 


C.3)  Model  of  Micro  turbines 

Two  types  of  microturbine  models  were  developed,  namely  (a)  the  detailed  model  and  (b)  the 
simplified  model.  The  detailed  model  was  developed  to  study  the  real-time  system  interaction  (in 
millisecond  timeframe).  However,  the  simulation  time  to  run  the  detailed  model  was 
significantly  long.  To  reduce  the  simulation  time,  the  simplified  model  was  developed  to  enable 
the  study  of  the  ID  APS  energy  management  system  for  a  24-hour  period. 

C.3.1)  The  detailed  microturbine  model 

a)  Model  Description: 

According  to  [69],  a  micro  turbine  model  comprised  a  gas  combustion  turbine  engine  integrated 
with  an  electrical  pennanent  magnet  synchronous  generator  (PMSG)  that  produced  electric 
power  while  operating  at  a  high  speed.  As  shown  in  Fig.  4,  the  inputs  of  the  gas  turbine  sections 
were  per  unit  turbine  speed  (output  from  PMSG)  and  a  rotor  reference  speed.  The  gas  turbine 
compared  the  rotor  speed  and  the  reference  speed  and  adjusted  its  outputs  —  per  unit  fuel  flow 
and  per  unit  turbine  torque  —  to  keep  the  rotor  speed  of  the  PMSM  at  the  reference  speed.  In 
addition  to  the  per  unit  fuel  demand  and  the  per  unit  turbine  torque,  the  other  output  of  the  gas 
turbine  section  was  the  exhaust  temperature  (F). 
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Fig.  4.  Microturbine  model  (Gas  turbine  subsystem  with  PMSG) 

As  shown,  torque  from  the  turbine  in  per  unit  was  converted  to  Newton-meter  (Nm),  and  then 
was  used  as  an  input  to  drive  the  PMSG.  Outputs  from  PMSG  were  three-phase  power  and  rotor 
speed  (rad/s).  The  rotor  speed  (rad/s)  was  converted  to  a  per  unit  quantity  and  was  fed  back  to 
the  turbine  section. 

b)  Mathematical  Formulae 

b.l)  The  gas  turbine  section: 

The  gas  turbine  section  of  this  model  was  built  according  to  the  mathematical  representations  of 
heavy-gas  turbines  presented  in  [70].  The  block  diagram  of  the  gas  turbine  section  is  shown  in 
Fig.  5.  The  validity  of  the  model  is  limited  to  simple  cycle  single-shaft,  generator  drive  gas 
turbine  only. 


Turbine  torque  (pu) 

3.  Compressor-turbine 


Fig.  5.  Gas  turbine  subsystems 
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The  dynamic  model  for  a  combustion  gas  turbine  consisted  of  the  following  subsystems:  (1) 
speed  control;  (2)  fuel  control;  (3)  compressor-turbine;  and  (4)  temperature  control.  All 
parameters  and  constants  used  in  the  gas  turbine  model  were  extracted  from  [71]. 

Speed  control  section:  The  speed  control  section  compared  the  current  rotor  speed  (from  PMSG) 
with  the  reference  speed  in  per  unit  quantities.  The  dynamic  response  of  the  speed  control  section 
was  modified  by  adjusting  parameters  in  the  transfer  function.  The  output  of  the  speed  control 
section  was  one  of  two  inputs  to  the  least  value  gate  (LVG).  The  other  input  was  from  the 
temperature  control  section.  The  lowest  value  of  the  two  inputs  was  sent  to  the  saturation  block 
(for  a  typical  gas  the  minimum  and  maximum  limits  are  -0.1  and  1.5,  respectively)  and 
forwarded  to  the  fuel  control  system. 

Fuel  con  trol  section  :  The  output  from  LVG  represented  the  least  amount  of  per  unit  fuel  needed 
for  a  particular  operating  point  that  corresponds  directly  to  the  per  unit  mechanical  power  on 
turbine  base  in  steady  state.  For  example,  if  mechanical  power  was  0.7  pu  then  the  steady  state 
value  for  fuel  demand  signal  was  also  0.7  pu.  To  obtain  the  final  fuel  demand  signal,  the  per  unit 
fuel  demand  from  the  LVG  was  multiplied  by  the  per  unit  turbine  speed.  The  output  was  then 
scaled  with  gain  K  =  0.77,  then  offset  by  the  constant  0.23,  which  represented  the  minimum 
amount  of  fuel  flow  at  no  load  rated  speed.  The  fuel  control  system  section  of  the  diagram  above 
could  be  interpreted  as  follows: 

•  When  no  more  fuel  required  (fuel  control  =  0),  fuel  demand  signal  =  0.23. 

•  When  fuel  control  =  1,  turbine  speed  =  0.5,  fuel  demand  signal  =  0.23+0.5*0.77  =  0.615. 

•  When  fuel  control  =  1,  turbine  speed  =  1.0,  fuel  demand  signal  =  0.23+1.0*0.77  =  1.0. 

The  fuel  control  section  also  consisted  of  the  fuel  valve  positioner  and  fuel  actuator.  The  transfer 
functions  associated  with  valve  positioner  and  fuel  actuator  generated  some  time  delays,  which 
in  this  case  was  modeled  at  0.05  and  0.04  seconds,  respectively. 

Compressor-turbine  section:  The  fuel  combustion  in  the  combustor  resulted  in  turbine  torque 
and  in  exhaust  temperature.  Inputs  of  the  compressor-turbine  section  were  the  per  unit  turbine 
speed  (N)  and  per  unit  fuel  flow  required  (Wf).  Outputs  of  the  compressor- turbine  section  were 
per  unit  torque  and  exhaust  temperature  (in  F),  given  as  inputs.  Torque  and  exhaust  temperature 
related  to  fuel  flow  (Wf)  and  turbine  speed  (N)  linearly  according  to  the  following  relationships. 

Torque  =  KHHV  (Wf -  0.23)  +  0.5  ( 1  -  N)  (Nm)  (Eq.  1 5) 

Exhaust  Temp  (TX)  =  Tref  -  700  (1  -  Wf)  +  550  (1  -  N)  (F)  (Eq.  16) 


Where: 

Wf  :  per  unit  fuel  flow 

N  :  per  unit  rotor  speed 

KHHV  :  coefficient  that  depends  on  the  higher  heating  value  of  the  gas 
stream  in  the  combustion  chamber  (KHHV  =1.2  for  a  typical 
gas) 

Tref  :  Reference  temperature  (950°F) 
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In  the  compressor-turbine  sections,  time  delay  associated  with  the  combustion  reaction,  time 
delay  associated  with  the  compressor  discharge  volume,  and  a  transport  delay  associated  with  the 
transportation  of  gas  from  the  combustion  system  through  the  turbine  were  set  at  0.01  second, 

0.2  second  and  0.04  second,  respectively. 

Temperature  control  section:  The  temperature  control  section  was  designed  to  limit  output 
power  of  the  gas  turbine  if  the  exhaust  temperature  of  the  turbine  was  higher  than  a  pre¬ 
determined  reference  temperature.  In  practice,  the  exhaust  temperature  was  measured  using  a 
series  of  thermocouples  incorporating  radiation  shields,  both  of  which  were  modeled  by  using 
transfer  functions  that  resulted  in  very  high  time  delays.  The  output  of  the  thennocouple  was 
compared  with  a  reference  temperature  (950F  in  this  case).  If  the  reference  temperature  remained 
higher  than  the  thermocouple  temperature,  which  was  usually  the  case,  this  would  permit  the 
dominance  of  speed  control  through  the  LVG.  When  the  thermocouple  output  exceeded  the 
reference  temperature,  the  difference  became  negative  and  the  temperature  control  output  would 
be  lower  than  the  speed  control  output.  This  allowed  the  temperature  control  signal  passing 
through  the  LVG  and  limits  the  turbine’s  output. 

b.2)  Permanent  Magnet  Synchronous  Generator 

A  pennanent  magnet  synchronous  generator  is  a  synchronous  generator,  in  which  the  DC  field 
winding  of  the  rotor  is  replaced  by  a  pennanent  magnet. 

The  voltage  rating  of  the  machine:  could  be  changed  by  varying  the  amount  of  flux  induced  by 
magnets.  The  relationship  between  flux  and  generate  voltage  of  the  machine  used  was: 

Ea  =  Km  ■  co(radls)  (Eq.  17) 


Where  Km  =  motor  torque  constant  (Nm/A),  which  is  a  function  of  flux  induced  by  magnet 
Ea  =  stator  voltage  (Vrms/phase) 
co  =  rotor  speed  (rad/s) 

K"'=^k  P  (P  (Eq'  18) 


p  =  number  of  poles 
cp  =  flux  (Wb) 

For  example,  if  we  would  like  to  design  a  PMSG  with  the  voltage  rating  of  480Vrms  (line-line) 
and  the  rotor  speed  of  70,000  rpm,  the  flux  per  phase  would  be  0.00891 1  Wb/phase. 

3 

Ea  =  —j=-  p-  cp-  a>(rad/s ) 

V2 

3  2  Jt 

=  -j=  •  2  •  0.00891 1  •  70000  x  — 

V2  60 

=  211.1  VrJ  phase  or  480  Vrmsil_l} 

Power  rating  of  the  machine:  could  be  changed  by  varying  the  resistance  and  reactance  of  the 
machine.  We  used  typical  resistance  and  reactance  values  for  a  synchronous  machine  of  0.01  per 
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unit  and  0.15  per  unit,  respectively.  For  example,  if  we  would  like  to  design  a  PMSG  with  the 
power  rating  of  30kW  and  the  voltage  rating  of  480Vrms  (line-line), 


Vbase  = 

Phase  = 

Ibase  = 

480/sqrt(3)  Vrms  (line-neutral) 

30kW 

36.1  A 

Zbase  = 

Vbase/Ibase 

=  7.68  Ohm 

Ra 

0.01x7.68 

=  0.0768  Ohm 

Xa 

0.15x7.68 

=  1.152  Ohm 

@  70,000rpm 

— ►  La  =  7.858x10-5  Henry 

C.3.2)  The  simplified  microturbine  model 

a)  Model  Description: 

A  simplified  micro  turbine  model  was  developed  following  the  characteristics  of  a  30kW 
capstone  microturbine,  in  terms  of  its  efficiency  and  the  rise  time.  As  shown  in  Fig.  6,  the 
model’s  input  was  the  microturbine’s  power  set  point  that  sets  the  amount  of  power  required 
from  the  microturbine.  The  output  of  the  model  was  its  power  output  (kW)  and  the  amount  of 
fuel  used  (BTU  or  cu  ft). 


Signs)  1 


-►  Agent  Ctr 


Set  point 


Pout 
Fuel  used 
ACout 


30kW  Microturbine 


£11 

Power  Meas. 

|-KlP2D 

Fuel  Meas. 

— VA— ° — 


R 


Fig.  6.  Microturbine  model  (Input  is  power  set  point;  outputs  are  power  output  and  fuel  and 

power  measurements) 

b)  Mathematical  Formulae  (Microturbine’s  Efficiency  and  its  Behavior  during  Load  Changes) 

Efficiency:  In  general,  the  efficiency  of  a  microturbine  decreases  when  it  is  partially  loaded.  At 
full  load,  the  efficiency  of  a  typical  30kW  microturbine  was  approximately  at  28%.  At  half  load, 
the  microturbine  efficiency  decreased  to  approximately  23%.  According  to  the  characteristic  of  a 
30-kW  microturbine  [72],  Fig.  7  illustrates  the  microturbine  part  load  power  performance. 
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Percent  Load  (%) 

Fig.  7.  Part  load  performance  (black)  of  a  30kW  microturbine  and  its  piecewise  linear  curve  fit 

(green)  (Source:  Energy  Nexus  Group) 

To  model  this  relationship  in  Matlab/Simulink,  the  piecewise  linear  curve  was  built  to  fit  the  part 
load  performance  curve.  This  was  also  shown  in  Fig.  7  as  a  solid  green  piecewise  linear  line. 

The  curve  fit  equations  used  were  as  follows: 

When  the  percent  load  (Li)  is  less  than  or  equal  to  15%,  r|i  =  (2*Li  +  0.35)  x  r|MT 

When  the  percent  load  (Li)  is  between  15%  and  40%,  rji  =  (0.8*Li  +  0.53)  x  r|MT 

When  the  percent  load  (Li)  is  between  40%  and  100%,  r|i  =  (0.25  *Li  +  0.75)  x  r|MT 

Where:  r|I  =  microturbine  efficiency  in  hour  i 

Li  =  load  in  percentage  of  microturbine  capacity  in  hour  i 
qMT  =  microturbine  efficiency  at  full  load  of  28% 

Rise  time:  Fig.  8  illustrates  the  characteristic  of  a  30kW  capstone  microturbine  when  responding 
to  the  power  demand.  Fig.  8  indicates  that  it  took  about  40  seconds  for  a  microturbine  to 
increase  the  power  output  from  zero  to  full  load  after  receiving  the  power  demand  signal.  This 
hot  start  characteristic  was  imitated  in  the  simplified  microturbine  model  using  a  signal 
generator. 
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Fig.  8.  Capstone  microturbine  test  -  cold  start  (Source:  CERT  [73]) 


Start-up  time:  The  start-up  characteristic  of  the  microturbine  is  also  illustrated  in  Fig.  9,  which 
implies  that  it  took  about  3  minutes  to  increase  the  turbine’s  power  output  to  full  load  from  its 
OFF  state.  Similar  to  the  cold  start  characteristic,  this  hot  start  characteristic  was  imitated  in  the 
simplified  microturbine  model  using  a  signal  generator. 


Fig.  9.  Capstone  microturbine  test  -  hot  start  (Source:  CERT  [73]) 


C.4)  Model  of  Fuel  Cells 


19 


The  fuel  cell  model  available  in  Matlab/Simulink  was  used  in  the  simulation.  The  mathematical 
fonnulae  for  the  model  development  were  extensively  discussed  in  Matlab  2009  Documentations 
and  would  not  be  repeated  here. 

C.5)  Model  of  Battery 

The  battery  model  available  in  Matlab/Simulink  was  used  in  the  simulation.  Slight  modifications 
were  made  internally  to  the  model  to  ensure  that  the  model  could  run  at  the  required  time  step 
(either  one  minute,  one  second  or  less).  The  mathematical  formulae  for  the  model  development 
were  extensively  discussed  in  Matlab  2009  Documentations  and  would  not  be  repeated  here. 


D)  Load  Models 

This  study  focused  on  household  loads  and  their  usage  characteristics.  Hourly  residential  load 
curves  of  an  average  household  used  in  this  study  were  extracted  from  the  RELOAD  database 
[74],  which  was  used  by  the  Electricity  Module  of  the  National  Energy  Modeling  System 
(NEMS).  The  hourly  residential  load  curve  data  were  available  for  twelve  months  (January  to 
December),  three  day  types  (typical  weekday,  typical  weekend  and  typical  peak  day)  and  nine 
load  types  (space  cooling,  space  heating,  water  heating,  cooking,  cloth  drying,  refrigeration, 
freezing,  lighting  and  others).  As  the  load  curves  in  the  RELOAD  database  represented  hourly 
fractions  of  the  yearly  load,  the  load  curves  were  scaled  up  by  the  annual  household  consumption 
and  divided  by  the  number  of  hours  in  a  year,  which  is  8760.  Therefore,  the  adjustment  made  to 
the  hourly  RELOAD  residential  load  curves  for  each  load  type  could  be  represented  by  Eq.  19. 

L  =  /  x  (Eq.  19) 

hou  8760 

where:  Lhour  =  Average  hourly  load  (kWh/h) 

f  =  Hourly  fraction  of  yearly  load 

Lannual  =  Average  annual  household  load 


4.2  Task  2:  Development  of  Local  Control  Algorithms  to  Control  ID  APS  Device  Models 

A)  Task  Description 

Task  2  involved  the  development  of  local  control  algorithms  to  control  DER  units,  and  loads. 
Particular  focus  was  given  to  the  DER  units  that  require  inverter  interface  with  the  utility  grid. 

B)  DER  Control 

Local  Control  during  the  Grid  Connected  Mode: 

During  the  grid  connected  mode,  voltage  magnitude  and  phase  at  the  interconnection  point  were 
the  two  control  parameters.  Voltage  magnitude  and  phase  from  DER  must  be  compatible  with 
those  at  the  interconnection  point  to  ensure  smooth  interconnection.  The  local  control  circuit  is 
illustrated  below  that  provides  the  grid  interface  between  a  DER  device  and  the  utility. 


Distributed! Grid  interface:  Inverter] 
generator  |  PWM  controller,  filter  j 


◄ - ► 


◄ 


Loads  & 
measurements 


Isolation  transformer 
&  the  grid 


► 
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◄ - ► 


Subsystems 


Fig.  10.  Test  circuit  for  the  grid  connected  operation 

As  shown,  the  grid  interface  section  comprises  an  inverter,  a  PWM  controller  and  a  filter.  The 
PWM  controller  is  responsible  for  calculating  voltage  magnitude  and  phase  based  on  the 
measurement  at  the  interconnection  point  (bus  B2  in  Fig.  10).  The  outputs  from  the  PWM 
controller  were  used  to  control  the  output  of  the  inverter  so  that  the  output  voltage  from  DER 
(bus  B 1  in  Fig.  10)  was  compatible  with  that  at  the  interconnection  point  (bus  B2)  at  all  time. 

The  PWM  controller  comprises  a  phase-lock  loop  (PLL)  unit  and  voltage  and  current  regulators 
(PI  controllers).  See  Fig.  11. 


Freq 

-H 

+1 

fabc  (pu)  wt 

Sir\_Cos 

* 

dqO 

►O 


*co 

phi(deg) 


Phase  Lead  2.4  deg 
compensating  phase  shift 
of  anti-aliasing  filters  at  60  Hz 


Fig.  1 1 .  PWM  controller 

The  following  regulator  gains  were  used  for  the  PLL  unit  and  the  voltage  and  current  regulators. 
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Proportional  gain 
(Kp) 

Integral  gain 

(Ki) 

Derivative  gain 
(Kd) 

PLL  block 

60 

1400 

0 

DC  voltage  regulator  block 

0.01 

3 

0 

Current  regulator  block 

0.05  (0.5) 

2(5) 

0 

Local  Control  during  the  Islanded  Mode: 

During  the  islanded  mode,  the  configuration  of  the  DER  control  unit  used  was  the  same  as  that 
during  the  grid  connected  mode.  In  this  case,  the  PWM  controller  would  monitor  the  voltage 
magnitude  and  phase  at  the  interconnection  point,  and  ensure  that  the  DER  could  connect  with 
the  rest  of  the  system.  For  a  system  with  a  single  inverter-based  DER,  the  PWM  setting  would  be 
set  to  follow  the  reference  voltage  at  1  per  unit  and  the  frequency  of  60  Hz  at  all  time. 

C)  Load  Control 

Load  control  strategies  were  designed  based  on  load  types.  In  this  study,  household  loads  were 
classified  into  three  groups,  namely  critical,  interruptible  and  deferrable  loads. 

•  Critical  loads  were  defined  as  loads  that  should  not  be  shut  down  or  shed  even  during 
emergency  condition. 

•  Deferrable  loads  were  referred  to  as  the  loads  that  could  be  shifted  to  operate  during  the 
later  hours  during  shortages  or  high  electricity  price  periods.  These  loads  may  include 
water  heating,  clothes  washing,  clothes  drying  and  other  plug-in  loads. 

•  Interruptible  loads  were  defined  as  the  loads  that  could  be  curtailed  during  emergencies 
or  during  the  high  electricity  price  periods.  These  loads  may  include  space  cooling  and 
optional  lighting  loads. 

The  load  control  strategy  was  developed  as  a  generalized  load  control  model  for  each  individual 
load  as  shown  in  Fig.  12. 


%  Compensation 
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Fig.  12.  Load  control  model  for  an  individual  load  developed  in  Matlab/Simulink 


As  shown,  the  model  required  the  followings  as  inputs:  (1)  price  signal  or  emergency  signal  from 
external  sources;  (2)  the  load  profile;  (3)  percentage  of  allowable  demand  reduction,  for 
example,  this  number  is  0%  for  critical  loads  and  is  100%  for  interruptible  loads;  and  (4) 
percentage  of  compensation  after  deferring  the  load,  for  example,  this  number  is  0%  for  critical 
loads  and  is  100%  for  interruptible  loads. 


4.3  Task  3:  Development  of  an  ID  APS  Energy  Management  System 

A)  Task  Description 

This  task  involved  the  development  of  the  ID  APS  energy  management  system  based  on  a  multi¬ 
agent  technology.  This  section  defines  the  environment  in  which  the  ID  APS  agents  resided,  the 
ID  APS  agent  architecture  —  which  included  types,  roles  and  their  interaction  —  and  lastly  the 
ID  APS  agent  decision  criteria  during  nonnal  operation  and  outage  durations. 

B)  Defining  the  Agent  Environment 

In  an  ID  APS  microgrid,  the  agents  interact  with  the  environment  through  sensors  and  actuators. 
Sensors  provide  agents  with  the  current  status  of  the  microgrid,  while  actuators  activated, 
controlled,  or  inactivated  a  process  or  action.  In  the  system  under  study,  ID  APS  agents  require 
sensors  to  sense  the  following  information  from  various  devices. 


Table  1.  Infonnation  sensed  and  calculated  by  agents 


Device/system 

Infonnation  Sensed 

Information  calculated 

DERs 

Magnitude  and  phase  of  voltage  and 
current  outputs 

Power  output  (kW),  fuel  used  (gallons 
or  BTUs),  emission  produced  (lbs) 

Battery 

Voltage,  charge/discharge  current, 
state  of  charge 

Power  charged  or  power  discharged 
(kW) 

Loads 

Magnitude  and  phase  of  voltage  and 
current  consumed  at  the  load 

Power  consumed  (kW) 

Switches 

Status  ON/OFF 

N/A 

Utility  system 

Electricity  price  ($/kWh)  -  the  price 
can  vary  from  every  minute  to  every 
hour;  system  frequency  and  voltage 

Whether  or  not  the  system’s  voltage  and 
frequency  are  in  the  safe  operation  range 

Agents  issue  the  following  control  signals  to  control  each  device  in  the  ID  APS  microgrid. 


Table  2.  Control  signals  from  agents 


Device/system 

Control  signals  from  agents 

DERs 

Power  set  point 

Battery 

Charge/discharge  signal 

Loads 

Status  ON/OFF,  consumption  reduction 
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Switches 

Status  ON/OFF 

Utility  system 

Disconnect  the  microgrid  once  an  upstream  failure  is  detected. 

C)  Defining  the  Agent  Architecture  (Type,  Role,  Interaction) 


Defining  Types  of  Agents: 

In  this  study,  the  idea  behind  any  multi-agent  system  is  to  break  down  a  complex  problem 
handled  by  a  single  entity  into  smaller  simpler  problems  handled  by  several  entities.  This  is 
called  a  distributed  and  decentralized  system.  The  proposed  multi-agent  system  consisted  of  four 
agents,  namely  control  agent,  DER  agent,  user  agent  and  database  agent.  The  four  agents 
together  formed  the  multi-agent  system  that  performed  actions  to  achieve  the  goal  of  the  system, 
that  is,  to  perfonn  demand  management  during  normal  operating  condition,  and  increase  the 
system’s  resilience  by  performing  islanding  and  securing  critical  loads.  The  architecture  of  an 
ID  APS  multi-agent  system  was  developed  and  presented  in  Fig.  13. 


Fig.  13.  ID  APS  agent  architecture  design  and  their  interaction,  where  arrows  represent 
messaging  exchange  among  agents  via  the  Transmission  Control  Protocol/Intemet  Protocol 

(TCP/IP) 

The  proposed  multi- agent  system  was  developed  to  follow  the  IEEE’s  standard  on  Foundation 
for  Intelligent  Physical  Agents  (FIPA).  This  helped  ensure  interoperability  among  different 
systems  and  platforms  so  that  the  proposed  multi-agent  system  would  be  universally  accepted. 


Defining  Agents  ’  Roles  and  responsibilities: 

•  The  control  agent :  monitoring  the  health  of  the  upstream  utility  network.  Once 
contingency  situations  or  grid  failures  are  detected,  it  sends  control  signals  to  isolate  the 
microgrid  from  the  utility.  It  is  also  the  control  agent  who  detects  the  restoration  of  the 
upstream  grid  and  issues  the  resynchronization  signal. 

•  The  DER  agent :  storing  the  associated  DER  information,  as  well  as  monitoring  and 
controlling  DER  power  levels  and  its  connect/disconnect  status.  DER  information  to  be 
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stored  may  include  DER  identification  number,  type  (solar  cells,  microturbines,  fuel 
cells,  etc),  power  rating  (kW),  local  fuel  availability,  cost  function  or  price  at  which  users 
agree  to  sell,  as  well  as  DER  availability,  i.e.  planned  maintenance  schedule. 

•  The  user  agent :  a  customer  gateway  that  makes  features  of  an  ID  APS  microgrid 
accessible  to  users.  A  user  agent  also  monitors  voltage,  current,  active  and  reactive  power 
consumption  at  each  critical  and  non-critical  load.  A  user  agent  also  allows  users  to 
control  the  status  of  loads  based  on  priority  pre-defined  by  a  user. 

•  The  database  agent :  storing  system  information,  as  well  as  recording  the  messages  and 
data  shared  among  agents.  Database  agent  also  serves  as  a  data  access  point  for  other 
agents,  and  keeps  track  of  all  available  agents  and  their  capabilities. 


Agents  ’  Collaborative  Diagram: 

This  section  summarizes  the  agent’s  collaborative  diagram  and  its  knowledge  modeling, 
including  Facts  and  Ontology.  The  collaborative  diagram  defines  the  interaction  among  agents 
and  their  interactions  with  the  environment.  For  example,  it  defines  how  the  control  agent 
interacts  with  the  distribution  circuit  and  its  measurement;  how  the  DER  agent  interacts  with 
DER  devices  and  measurements;  and  how  the  user  agent  interacts  with  users  and  load 
measurements.  The  collaborative  diagram  also  helps  fonnulate  the  agent  roles  and 
responsibilities,  which  simplify  understanding  and  modeling  the  problem  at  hand. 

Fig.  14  represents  the  ID  APS  agent’s  collaborative  diagram,  which  illustrates  three  ID  APS 
agents  -  a  control  agent,  a  DER  agent  and  a  user  agent  -  and  their  interactions  with  each  other 
and  the  environment.  All  three  ID  APS  agents  communicate  with  the  ID  APS  database  agent, 
which  comprises  the  name  server  and  the  facilitator.  The  visualizer  receives  copies  of  all 
messages  exchanged  within  the  multi-agent  system  and  is  responsible  for  displaying  these 
messages  to  users.  In  this  study,  the  name  server,  the  facilitator  and  the  visualizer  are  collectively 
termed  utility  agents. 


Fig.  14.  Agent’s  collaborative  diagram 
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Numbers  1-22  represent  interactions  among  agents  or  between  the  agents  and  their  environment. 
These  associated  interactions  are  summarized  in  Table  3. 

Table  3.  Various  agents’  actions  in  the  collaborative  diagram 


Collaboration  Explanation 


1 

Registration 

Agents  notify  the  database  agent  of  their  presence 

2 

Location  and  type  query 

Request  for  the  network  location  of  agents  and  types 

3 

Location  and  type 
response 

Response  from  agents  for  their  locations  and  types 

4 

Ability  request 

Request  for  agents’  abilities 

5 

Ability  response 

Response  from  agents  about  their  abilities 

6 

Find  response 

Response  containing  a  list  of  agents  matching  the  desired  criteria 

7 

Find  request 

Request  to  obtain  a  list  of  agents  with  particular  abilities 

8 

Activity  notification 

Response  containing  all  messages  sent  among  agents  or  between 
agents  and  their  environment 

9 

Activity  request 

Request  to  forward  all  activities 

10 

UA  registration  request 

Request  from  user  agent,  to  control  agent  for  subscription 

11 

DERA  registration  req 

Request  from  DER  agent,  to  control  agent  for  subscription 

12 

CA  registration  response 

Ack.  from  control  agent  that  user  agent’s  ID  has  been  stored. 

13 

CA  registration  response 

Ack.  from  control  agent  that  DER  agent’s  ID  has  been  stored. 

14 

CA  information  retrieval 

Retrieval  of  information  from  the  main  grid  by  the  control  agent 

15 

CA  inform  signals 

Inform  signals  issued  by  the  control  agent  to  inform  all  registered 
subscribers  about  the  information  retrieved  from  the  main  grid 

16 

CA  control  signals 

Control  signals  issued  by  control  agent  to  open/close  main  circuit 
breaker 

17 

UA  information  retrieval 

Retrieval  of  loads’  power  requirement  from  external  circuit  by  user 
agent 

18 

UA  load  information 

Information  sent  by  user  agent  to  inform  DER  agent  on  load  power 
req 

19 

DERA  info  retrieval 

Retrieval  of  information  from  DERs  by  the  DER  agent 

20 

DER  output  information 

Information  sent  by  DER  agent  on  amount  of  power  sent  to  user 
agent 

21 

UA  control  signals 

Control  signals  issued  by  the  user  agent  to  turn  the  load  ON  or  OFF 

22 

DERA  control  signal 

Control  signals  issued  by  DER  agent  to  control  power  production 
from  DER 

Defining  Agents  ’  Knowledge  Modeling:  Facts  and  Ontology: 

In  knowledge  modeling,  Facts  and  Ontology  were  defined  for  each  agent’s  application.  Facts 
represent  statements  that  an  agent  believes  to  be  true,  either  about  itself  or  its  external 
environment.  Ontology  is  the  vocabulary  used  by  agents  during  communication.  Ontology  is 
created  for  the  multi-agent  system  based  on  the  application  design  process  mentioned  in  the 
previous  section.  The  application  consists  of  both  abstract  and  physical  concepts.  For  instance, 
islanding  mode  information  is  abstract,  while  Energy  is  physical  concept.  Table  4  lists  Facts 
defined  in  the  proposed  multi-agent  system,  together  with  their  attributes  and  default  values. 

Table  4.  Facts  in  the  ID  APS  multi-Agent  system 
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Facts 

Additional  Attributes 

Default  Value 

islandmode 

is  island:  Boolean 
main  cb  status:  Boolean 
id:  String 

outageDuration:  String 
outageHour:  String 

is  island:  false 
main  cb  status  :true 
id:  n/a 

outageDuration:  n/a 
outageHour:  n/a 

criticalCB 

clcb  status:  Boolean 

clcb  status:  true 

noncriticalCB 

nclcb  status:  Boolean 

nclcb  status:  true 

DERspecs 

sendDERinfo:  Boolean 

DER  Cost:  String 
UserAgentNameSpecs:  String 

sendDERinfo:  false 

DER  Cost:  “0” 
UserAgentNameSpecs:  n/a 

ua  DERcmd 

der_required_power:  String 
send_der_cmd:  Boolean 
UserAgentName:  String 

der  require  power:  0.0 
send_der_cmd:  false 

UserAgentName:  n/a 

re  uaDERcmd 

required_power:  String 
UserAgentName:  String 

require  power:  0.0 

UserAgentName:  n/a 

dgDERcmd 

der_produced_pwr:  String 
send_der_cmd:  Boolean 
userAgents:  String 

der  produced  pwr:  0.0 
send_der_cmd:  false 
userAgents:  n/a 

redgDERcmd 

producedjpwr:  String 
sendDERCntrlCmd:  String 

producedjpwr:  0.0 
sendDERCntrlCmd:  false 

agentsName 

name:  String 

name:  n/a 

agentRegistered 

name:  String 

name:  n/a 

inIDAPS 

name:  String 
isInIDAPS:  Boolean 

name:  n/a 
_is!n!DAPS:  false 

More  detailed  information  on  this  topic  can  be  found  in  [75,  76]. 


D)  Agent  Decision  Criteria 

During  Normal  Operation: 

During  normal  operations,  the  overall  objective  of  multi-agent  collaboration  was  designed  to 
achieve  look-ahead  heterogeneous  dispatch  by  controlling  both  loads  and  generator.  The  loads 
and  generators  were  controlled  separately. 

•  Generator  control:  DERs  were  dispatched  according  to  their  merit  order  in  responding  to 
the  fluctuation  in  the  demand  for  electricity.  The  entire  amount  of  energy  generated  from 
renewable  energy  sources  such  as  solar  cells  or  wind  turbines  were  injected  into  the 
distribution  circuit.  Mixed-integer  linear  programming  was  used  to  solve  this  scheduling 
problem.  The  decision  variables,  the  objective  function  and  the  constraints  are  described 
below. 

Decision  variables: 

DER  output  (Xij)  to  be  dispatch  at  each  hour  of  the  day 
Objective  function: 

Minimize  the  operating  cost  at  each  hour  of  the  day. 
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Minimize 


2(<vX) 

1=1 

where 

Xi :  Decision  variables,  power  set  point  (kW)  of  DER  i 

ay.  Operating  cost  ($/kWh)  of  DER  i 

Constraints: 

1.  All  loads  must  be  served  on  demand,  except  for  the  deferrable  ones. 

n 

For  each  hour,  ^(W)  >  L  L  is  load  requirement  (kW) 

/=  1 

2.  The  cumulative  fuel  used  must  still  be  less  than  the  available  fuel  as  indicated  by  fuel 
sensors. 

For  each  hour,  0  <  Xt  <  C,  ■  Fi 

Where: 

Ci:  Capacity  of  DER  i  (kW) 

Fi:  Fuel  availability  indicator  of  DER  i,  if  fuel  is  available  Fi  =  1,  otherwise  Fi  =  0 


3.  The  cumulative  emission  produced  must  still  be  less  than  the  allowable  emission  limit. 
For  each  hour,  0  <  Xt  <  C,  •  E{ 

Where: 

Ci:  Capacity  of  DER  i  (kW) 

Ei:  Emission  limit  indicator,  if  the  total  emission  is  still  within  the  available 

range  Ei  =  1 ,  otherwise  Ei  =  0 

•  Load  control:  during  the  grid-connected  mode,  the  load  control  strategy  was  set  to  be 
responsive  to  the  real-time  electricity  price  during  a  non-emergency  condition.  Load 
control  strategies  were  discussed  in  Section  4.2. 


During  Outage  Conditions: 

During  upstream  outage  conditions,  the  ID  APS  microgrid  was  designed  to  disconnect  itself  from 
the  utility  grid  and  secure  critical  loads  according  to  the  given  prioritized  list.  Some  non-critical 
loads  could  also  be  secured  depending  upon  the  user’s  preference.  During  each  simulation,  the 
first  step  was  to  identify  critical  loads  and  their  priority.  That  is,  each  homeowner  was  asked  to 
identify  their  critical  loads  to  be  secured  during  an  event  of  an  upstream  outage.  Priority  of  these 
critical  loads  was  also  identified.  If  available  generation  was  not  sufficient  to  secure  all  critical 
loads,  the  least  important  critical  loads  were  shed.  In  the  simulation,  a  homeowner  could  also 
specify  the  time  of  day  when  he/she  would  like  to  operate  their  non-critical  loads  during 
emergencies.  As  long  as  there  was  excess  capacity  from  emergency  generators,  and  expenses 
required  to  run  these  loads  were  within  a  reasonable  range,  selected  non-critical  loads  would  be 
secured  during  the  chosen  duration. 
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4.4  Task  4:  Development  IP-based  Communication  Interfaces  for  Connecting  between  the 
ID  APS  Physical-Cyber  Components 

A)  Task  Description 

The  objective  of  this  task  was  to  demonstrate  how  a  connection  between  the  ID  APS  hardware 
elements  and  the  ID  APS  multi-agent  decision  support  systems  could  be  established  using  TCP/IP 
Client/Server  socket  communications. 

B)  TCP/IP  Connection 

In  the  simulated  environment,  the  multi-agent  system  connected  to  the  microgrid  in  the 
MATLAB/Simulink  environment  over  a  TCP/IP  connection.  This  allowed  the  microgrid  to  be 
controlled  remotely  by  the  multi-agent  system  from  any  location.  A  third  party  TCP  server 
implementable  in  MATLAB  Simulink  was  used  to  establish  the  TCP  connectivity.  Socket 
programming  was  carried  out  in  agents’  external  java  classes.  The  third  party  TCP  server  was 
capable  of  managing  62  outputs  and  unlimited  number  of  inputs.  The  server  allowed  only  a 
single  TCP  connection  at  a  time  and  followed  a  specific  format  for  the  input/output  messages. 
This  limited  our  control  over  the  connection  between  the  multi-agent  system  and  the  microgrid 
because  each  agent  (a  user  agent,  a  DER  agent  and  a  control)  required  a  separate  TCP  connection 
to  the  microgrid.  To  handle  the  situation  a  middle  server  was  developed. 

C)  Middle  Server 

The  middle  server  was  developed  to  allow  multiple  TCP  connections  to  the  microgrid  and 
flexible  control  over  the  message  exchange  between  the  microgrid  and  the  multi-agent  system.  It 
acted  as  a  server  for  the  multi-agent  system  and  a  client  for  the  microgrid’s  TCP/IP  server.  It 
allowed  multiple  input  TCP  connections  from  the  multi- agent  system  and  established  a  single 
TCP  connection  with  the  microgrid’s  TCP/IP  Server.  Any  message  sent  from  the  multi-agent 
system  was  processed  in  the  middle  server  and  sent  out  to  the  microgrid’s  TCP/IP  server  in  the 
appropriate  format.  Similarly,  messages  sent  by  the  microgrid’s  TCP/IP  server  were  received  by 
the  middle  server  and  sent  out  to  the  appropriate  receiver  in  the  multi-agent  system.  The  figure 
below  depicted  this  scenario. 
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Fig.  15.  Middle  server  implementation 


4.5  Task  5:  Simulation  and  Evaluation  of  the  Microgrid  in  both  Parallel  and  Islanded 
Operations 

A)  Task  Description 

In  order  to  test  this  idea,  we  simulated  and  analyzed  the  proposed  ID  APS  microgrid  concept 
using  realistic  data  from  Virginia  Tech  Electric  Service  (VTES).  The  Virginia  Tech  local 
distribution  network  was  an  ideal  example  of  an  interconnected  microgrid  in  a  small  university 
town. 

B)  Assumptions 

The  following  assumptions  were  made  during  the  simulation: 

•  The  electric  power  infrastructure  was  available  to  support  load  and  generator  control. 

•  The  communication  infrastructure  was  available  to  support  the  communications  among 
various  physical  devices  in  the  microgrid. 

•  Sufficient  DER  capacity  was  available  to  secure  critical  loads. 

•  The  IP-addressable  interface  was  available  to  support  connections  among  all  ID  APS 
physical  devices  and  the  agents. 

C)  Scenario  Description 

A  section  of  a  distribution  circuit  from  VTES  was  used  as  a  basis  for  this  simulation  study.  The 
ID  APS  microgrid  was  simulated  in  the  following  scenarios  to  demonstrate  its  ability  to  operate 
in  parallel  and  islanded  operations. 

Scenario  1 :  Demand  management  during  nonnal  operation 

The  objective  of  this  scenario  was  to  showcase,  during  the  non-emergency  operation,  that  the 
ID  APS  energy  management  system  could  perform  load  control  and  the  level  of  load  control 
depended  on  hourly  electricity  prices  and  load  types.  The  results  shown  in  Section  5  —  the  results 
and  discussion  section  —  illustrated  that  the  ID  APS  demand  management  features  would  reduce 
the  peak  demand,  resulting  in  the  reduction  in  demand  charges  from  a  local  electric  utility. 

Scenario  2:  Islanding  and  securing  critical  loads  during  an  emergency 

The  objective  of  this  scenario  was  to  illustrate  the  ability  of  multiple  ID  APS  agents  working  in 
collaboration  to  operate  the  microgrid  in  the  islanded  operation.  The  ID  APS  microgrid  was 
simulated  for  24  hours,  assuming  that  the  outage  occurred  at  5:00PM  for  three  hours.  The 
simulation  also  focused  on  how  the  ID  APS  agents  detected  the  upstream  fault,  as  well  as  how 
voltage  and  frequency  changed  during  islanding  and  re-synchronizing  of  the  microgrid  after  the 
upstream  fault  was  cleared. 

Scenario  3:  Serving  non-critical  loads  during  an  emergency 

The  objective  of  this  scenario  was  to  showcase  that  the  excess  capacity  from  the  DER  unit  after 
serving  critical  loads  could  be  shared  among  non-critical  loads  in  the  community.  As  a  result,  a 
DER  unit  can  be  fully  utilized  in  any  given  emergency  event. 
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5.  Results  &  Discussions 


All  tasks  were  completed  as  planned.  The  project  accomplishments  were  summarized  below. 


Tasks 

Tasks  accomplished 

1 .  Development  of  ID  APS 
device  models 

Developed/modified  simulation  models  for  PV  generators, 
wind  turbines,  microturbines,  fuel  cell  and  battery  storage 
in  Matlab  for  used  in  the  simulation;  Developed  load 
models. 

2.  Development  of  local 
control  algorithms  to  control 

ID  APS  device  models 

Developed  local  control  for  DER  devices  and  load  control 
strategies. 

3.  Development  of  an  ID  APS 
energy  management  system 

Defined  the  agent  environment,  the  agent  architecture,  and 
the  agent  decision  criteria 

4.  Development  of  IP-based 
communication  interface 

Established  a  TCP/IP  connection  between  the  ID  APS 
hardware  elements  and  the  ID  APS  decision  support  system 
in  a  simulated  environment 

5.  Simulation  and  evaluation 
of  the  microgrid  in  both 
parallel  and  islanded 
operations 

Simulated  the  ID  APS  microgrid  in  parallel  operation; 
Simulated  the  ID  APS  microgrid  in  islanded  operation; 
Evaluated  the  results  and  submitted  the  final  report 
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5.1  Task  1:  Development  of  ID  APS  Device  Models 


A)  PV  Generators  -  Sample  Model  Outputs 

The  model  of  a  PV  generator  was  developed  in  Matlab/Simulink.  Fig.  16  illustrates  the  model 
mask  captured  from  Matlab,  together  with  the  description  of  the  model’s  inputs/outputs. 


PV  Model 


As  shown  in  Fig.  16,  inputs  required  include: 

-  Open-circuit  voltage  (Voc), 

-  Number  of  cells  connected  in  series  (Ns), 

-  Short-circuit  current  (Isc), 

-  Number  of  cells  connected  in  parallel  (Np), 

-  Maximum  module  power  (Pmax), 

-  Temperature  coefficient  of  open-circuit 
voltage,  and 

-  Temperature  coefficient  of  short-circuit 
current. 

These  parameters  can  be  easily  obtained  from  a 
manufacturer’s  datasheet.  Additional  inputs  to 
the  model  are  cell  temperature  and  solar 
irradiation,  which  can  be  obtained  from 
historical  data. 


Fig.  16.  PV  model 

The  model  outputs  were  validated  by  comparing  the  I-V  curve  outputs  with  the  I-V  curves  given 
by  the  manufacturer.  As  shown  below,  the  I-V  curves  generated  by  the  developed  model  were 
very  similar  to  those  of  the  manufacturer,  using  the  parameters  of  a  BP-4175  175W  PV  module. 


Fig.  17.  The  I-V  curve  from  the  manufacturer  for  BP-4175  175W  PV  module  (left);  the  I-V 
curve  generated  by  the  developed  model  (right) 
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B)  Wind  Turbines  -  Sample  Model  Outputs 

As  discussed  in  Section  4.1,  a  wind  turbine  generator  was  modeled  as  a  wind  turbine  section 
connected  with  an  electrical  generator  section.  Fig.  18  illustrates  the  model  mask  captured  from 
Matlab.  Model  inputs  were  (1)  the  generator  speed  in  per  unit  of  the  nominal  speed  of  the 
generator,  and  (2)  the  wind  speed  in  m/s.  The  output  of  the  wind  turbine  model  was  the  three- 
phase  power  that  could  be  directly  connected  to  the  electrical  grid. 


Fig.  18.  Wind  turbine  model  =  wind  turbine  +  asynchronous  generator 

The  outputs  of  the  developed  wind  turbine  model  are  presented  in  Fig.  19.  Fig.  19  (top)  displays 
the  active  power  (yellow)  and  reactive  power  (pink)  outputs  of  the  developed  wind  turbine  model 
in  kW  and  kVAR,  respectively,  when  the  wind  speed  decreases  from  15m/sto  lOm/sin  the 
middle  of  the  simulation.  The  power  output  decreases  as  the  wind  speed  decreases.  The  wind 
turbine  absorbs  reactive  power  while  producing  electrical  power  output.  Fig.  19  (bottom) 
displays  the  three-phase  electrical  output  —  voltage  (yellow)  and  current  (pink)  —  in  per  unit. 
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Fig.  19.  Electrical  outputs  of  wind  turbine: 
Top:  Power  (yellow)  and  Reactive  Power  (pink); 
Bottom:  Voltage  (yellow)  and  Current  (pink) 
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C)  Microturbines  -  Sample  Model  Outputs 

As  discussed  in  Section  4.1,  two  microturbine  models  were  developed:  the  detailed  and  the 
simplified  models.  Fig.  20  illustrates  the  sample  outputs  of  the  detailed  microturbine  model,  e.g. 
the  line-to-line  voltage  of  a  30kW  microturbine  at  480V  line-to-line. 


Line-to-line  voltage  (Volt)  during  the  hot  start 
when  simulating  from  0  to  50  seconds 


Snap  shot  of  the  line-to-line  voltage  (Volt) 
around  t  =  50  second 


Fig.  20.  Sample  outputs  of  the  detailed  microturbine  model 
x-axis  is  the  simulation  time  in  second;  y-axis  is  the  voltage  magnitude; 


Fig.  21  illustrates  the  sample  outputs  of  the  simplified  microturbine  model,  e.g.  the  power  output 
of  a  30kW  microturbine  during  the  hot  start  (left)  and  the  cold  start  (right). 
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Fig.  21.  Sample  outputs  of  the  simplified  microturbine  model 
Characteristics  during  the  hot  start  (left);  Characteristics  during  the  cold  start  (right) 


D)  Fuel  Cells  and  Battery  Storage  Model 

Since  the  fuel  cell  and  the  battery  storage  models  available  in  Matlab/Simulink  were  used  in  the 
simulation,  the  model  outputs  could  be  found  in  Matlab  2009  reference  guide.  Therefore,  the 
sample  results  of  these  two  models  are  not  discussed  here. 
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E)  Load  Profiles  -  Sample  Model  Outputs 

Fig.  22  and  Fig.  23  illustrate  the  output  of  the  load  models  developed,  i.e.  load  profiles  of  a 
typical  home  in  summer  and  winter,  respectively. 
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Fig.  22.  Load  profile  of  a  typical  home  in  summer  (August) 
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Fig.  23.  Load  profile  of  a  typical  home  in  winter  (January) 
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5.2  Task  2:  Development  of  Local  Control  Algorithms  to  Control  ID  APS  Device  Models 
A.  Grid  Connected  Mode  -  Sample  Results 

As  shown  in  Fig.  24,  the  test  circuit  comprises  one  distributed  generator,  a  grid  interface  section 
-  an  inverter,  a  PWM  controller  and  an  RLC  filter,  loads  and  measurements,  as  well  as  an 
isolation  transformer  and  the  grid. 


Fig.  24.  Test  circuit  for  the  grid  connected  operation 

The  circuit  as  shown  was  simulated  for  0.2  seconds  to  illustrate  the  impact  of  the  decrease  in 
DER  output  on  the  voltage  signal  at  the  interconnection  point.  Fig.  25  displays  the  simulation 
results  when  the  power  output  from  a  DER  unit  decreases  from  50  kW  to  25  kW  at  t  =  0. 1 
seconds. 
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Fig.  25.  Simulation  result  illustrating  smooth  voltage  and  current  output  at  the  interconnection 
point  (bus  B2)  with  the  change  in  DER  power  output.  Where: 

Channel  1 :  voltage  and  current  of  phase  a  monitored  at  bus  B2  (pu) 

Channel  2:  voltage  of  phase  a  monitored  at  bus  B1  (volt) 

Channel  3:  power  command  (pink)  vs  power  output  (yellow)  measured  at  the  load  (kW) 

Channel  4:  DC  reference  voltage  (pink)  vs  DC  voltage  measured  at  the  capacitor  (yellow) 
(volt) 

The  results  illustrates  that  the  implemented  DER  control  strategy  results  in  smooth 
measurements  (voltage  &  current)  at  the  interconnection  point,  even  with  the  change  in  DER 
power  output. 

B.  Islanded  Operation  Mode  -  Sample  Results 

The  test  circuit  used  to  illustrate  the  local  control  model  during  the  islanded  mode  is  shown 
below.  This  is  the  similar  circuit  as  shown  in  Fig.  26  with  an  additional  fault  unit  at  the  utility 
side.  In  this  case,  it  is  assumed  that  the  upstream  fault  occurs  at  0.1  second.  After  the  fault 
occurs,  the  outputs  of  the  distributed  generator  were  controlled  to  follow  the  reference  voltage 
and  frequency. 
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Fig.  26.  Test  circuit  for  the  islanded  operation 


As  shown  in  Fig.  27,  the  control  model  implemented  allows  the  smooth  transition  from  the  grid 
connected  to  the  islanded  mode  of  operation  (at  t=0.1  second),  as  well  as  allowing  the  microgrid 
to  operate  autonomously  during  the  islanded  mode  (at  t  beyond  0. 1  second). 


Fig.  27.  Simulation  result  illustrating  smooth  voltage  and  current  output  at  the  interconnection 

point 


5.3  Task  3:  Development  of  an  ID  APS  Energy  Management  System 

Under  this  task,  the  multi-agent  system  for  managing  the  ID  APS  microgrid  was  successfully 
developed.  The  results  of  this  task  are  integrated  in  Section  5.5,  which  discusses  the  simulation 
results  when  the  agents  control  the  microgrid  in  both  parallel  and  islanded  operations. 
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5.4  Task  4:  Development  IP-based  Communication  Interfaces  for  Connecting  between  the 
ID  APS  Physical-Cyber  Components 

The  IP-based  communication  interface  to  connect  the  ID  APS  physical-cyber  systems  was 
successfully  developed.  The  results  of  this  task  are  presented  with  the  simulation  results  of  Task 
5  that  allows  agents  to  communicate  with  physical  elements  in  the  simulation  circuit  via  TCP/IP. 


5.5  Task  5:  Simulation  and  Evaluation  of  the  Microgrid  in  both  Parallel  and  Islanded 
Operations 

As  previously  discussed  in  Section  4.5,  the  ID  APS  microgrid  was  simulated  in  the  following 
three  scenarios  to  demonstrate  its  ability  to  operate  in  parallel  and  islanded  operations. 

•  Scenario  1 :  Demand  management  during  normal  operation 

•  Scenario  2:  Islanding  and  securing  critical  loads  during  an  emergency 

•  Scenario  3:  Serving  non-critical  loads  during  an  emergency 
The  simulation  results  are  discussed  below: 


A)  Case  study  1:  Perform  Demand  Management  during  Non-Emergency  Situations 

To  demonstrate  the  ability  of  the  ID  APS  agents  to  perform  load  control  in  response  to  the  price 
change,  we  used  electricity  tariffs  from  two  electric  utilities  in  the  Washington  DC  area: 
Baltimore  Gas  &  Electricity  (BG&E)  and  Dominion  Virginia  Power  (DOM). 

Demand  response  in  a  summer  month 

Fig.  28  (left)  illustrates  the  TOU  rates  of  both  BG&E  and  DOM  in  August  (the  summer  peak 
month).  The  prices  shown  here  represented  the  companies’  energy  and  delivery  charges  in  2009 
and  did  not  include  other  service  charges  and  taxes.  Fig.  28  (right)  illustrates  the  original 
household  load  shape  (gray  area)  and  the  load  shapes  after  implementing  demand  response 
strategies  (solid  blue  and  red  lines)  under  BG&E’s  and  DOM’s  TOU  rates. 


|  □  Load  BGE  Dominion  | 


Summer  Load  under  Different  TOU  Rates 


Fig.  28.  The  summer  TOU  rates  from  the  chosen  utilities  (left);  Summer  load  profile  (right)  - 
before  (gray  area)  and  after  the  demand  response  (blue  and  red  lines) 
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As  shown,  the  demand  response  strategies  implemented  by  ID  APS  agents  helped  lower  the 
system  peak  demand,  under  both  BG&E  and  DOM  rate  structures.  This  resulted  in  the  reduction 
in  a  portion  of  the  demand  charge  (S/kW)  imposed  by  a  utility.  With  the  ID  APS  demand 
management  strategy,  lower  system  peaks  were  expected  under  the  pricing  scheme  with  higher 
variation  in  peak  and  off-peak  rates.  In  this  case,  the  level  of  demand  response  would  be  higher 
under  DOM’s  TOU  rates,  as  compared  to  that  under  BG&E’s  TOU  rates. 


Demand  response  in  a  winter  month 


Fig.  29  (left)  illustrates  the  TOU  rates  of  both  BG&E  and  DOM  in  January  (the  winter  peak 
month).  The  prices  shown  here  represented  the  companies’  energy  and  delivery  charges  in  2009 
and  did  not  include  other  service  charges  and  taxes.  Fig.  29  (right)  illustrates  the  original 
household  load  shape  (gray  area)  and  the  load  shapes  after  implementing  demand  response 
strategies  (solid  lines)  under  BG&E’s  and  DOM’s  TOU  rates. 
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Fig.  29.  The  summer  TOU  rates  from  the  chosen  utilities  (left);  Summer  load  profiles  (right)  - 
before  (gray  area)  and  after  the  demand  response  (blue  and  red  lines) 

Similar  to  the  results  in  the  summer  month,  the  demand  response  strategies  implemented  by 
ID  APS  agents  could  help  lower  the  system  peak  demand,  under  both  BG&E  and  DOM  rate 
structures  in  the  winter  month  as  well.  More  detailed  description  of  ID  APS  control  algorithms 
during  non-emergencies  can  be  found  in  [77,  78]. 

B)  Case  study  2:  Island  the  Microgrid  once  a  Fault  is  Detected  and  Secure  Critical  Loads 
during  an  Emergency 

This  scenario  illustrated  the  ability  of  multiple  agents  to  work  in  collaboration  to  operate  the 
microgrid  in  the  islanded  mode.  The  ID  APS  microgrid  was  simulated  for  24  hours,  assuming  that 
the  outage  occurred  at  5:00PM  for  three  hours.  The  simulation  also  focused  on  how  ID  APS 
agents  detected  the  upstream  fault,  voltage  and  frequency  deviation  during  islanding  and  re¬ 
synchronizing  of  the  microgrid  after  the  upstream  fault  was  cleared. 

Fig.  30  illustrates  the  load  profiles  of  one  house  in  this  microgrid.  The  result  clearly  shows  that 
as  soon  as  an  outage  was  detected  the  multi-agent  system  performed  its  actions  to  secure  critical 
loads  prioritized  according  to  the  homeowner’s  choice.  In  this  case,  the  critical  loads  were 
secured  and  non-critical  loads  were  disconnected. 
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Fig.  30.  Load  profile  with  agents  securing  critical  loads  during  the  outage  between  5-8pm.  All 

non-critical  loads  are  shed. 

The  upstream  fault  was  applied  to  the  circuit  at  5  PM.  As  soon  as  upstream  fault  was  detected, 
the  ID  APS  agents  islanded  the  microgrid.  The  result  shown  in  Fig.  3 1  is  a  snapshot  of  the  event 
at  around  5  PM.  This  illustrates  the  fluctuation  in  voltage  and  current  at  the  interconnection  point 
of  the  ID  APS  microgrid  and  the  utility  system,  from  0. 1  seconds  before  the  fault  and  0.2  seconds 
after  the  fault.  This  helps  determine  how  quickly  the  multi-agent  system  can  perform  its  actions 
to  isolate  the  microgrid  while  maintaining  the  system  stability. 

pu 
l.o 
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-l.o 

5PM-0.1sec  5PM  5PM+0.1sec  5PM+0.2sec 

time  Voltage 

— Current 

Fig.  3 1 .  The  variation  of  60Hz  voltage  and  current  waveforms  measured  at  the  interconnection 
point  before  and  after  the  upstream  outage  applied  to  the  circuit. 

A.  During  the  Grid-Connected  Mode  (before  5PM) 

The  microgrid’s  voltage  and  frequency  were  controlled  such  that  they  followed  the  grid’s  voltage 
and  frequency,  roughly  at  1  per  unit  and  60  Hz,  respectively.  As  shown,  the  voltage  (yellow)  was 
always  kept  at  around  1  per  unit.  The  current  (pink)  indicated  the  current  produced  by  the  DER. 
Before  the  upstream  fault  occurred,  the  DER  unit  supplied  some  power  to  the  microgrid  at  about 
50%  of  its  rated  power. 

B.  During  the  Transition  (at  5PM) 

Once  the  control  agent  detected  the  fault  at  5PM,  the  control  agent  informed  the  user  agent  and 
the  DER  agent,  both  of  which  exchanged  information  and  determined  the  amount  of  loads  to  be 
shed  and  the  amount  of  power  to  be  produced  internally  in  order  to  stabilize  the  microgrid.  The 
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control  agent  sent  a  control  signal  to  isolate  the  microgrid  by  disconnecting  the  main  circuit 
breaker  as  soon  as  the  upstream  outage  was  detected.  Depending  on  the  pre-defined  load  priority 
set  by  a  user  and  the  available  internal  generation,  the  user  agent  disconnected  non-critical  loads 
and  certain  critical  loads.  In  this  case,  we  assumed  that  the  DER  could  fully  support  connected 
critical  loads  without  any  energy  limitation  during  the  outage.  In  other  words,  the  DER  could 
increase  its  output  instantaneously  with  the  help  of  battery  storage  once  the  microgrid  islanded. 
Thus,  the  user  agent  only  disconnected  the  non-critical  loads  and  left  the  critical  loads  intact.  All 
agent  actions  -  from  detecting  the  fault,  disconnecting  the  main  circuit  breaker,  disconnecting 
the  non-critical  loads  to  stabilizing  the  grid  -  were  accomplished  within  about  half  an  electrical 
cycle,  i.e.  less  than  0.008  second  for  a  60-Hz  system.  The  DER  fault  current  as  shown  was 
limited  by  the  inverter  rating  in  the  case  of  an  inverter-based  DER. 

C.  During  the  Islanded  Mode  (5PM-8PM): 

After  the  fault  occurred  at  5PM,  the  ID  APS  microgrid  was  disconnected  from  the  main  grid. 
When  the  microgrid  operated  in  the  island  mode,  the  grid  interface  unit  tried  to  control  the  DER 
output  voltage  at  1  per  unit  and  the  frequency  at  60Hz.  The  user  agent  and  the  DER  agent 
balanced  the  demand  and  supply  within  the  microgrid.  The  total  system  load  was  reduced  to  its 
critical  loads  and  the  DER  produced  almost  at  its  full  capacity  internally  to  supply  these  critical 
loads.  This  can  be  seen  in  Fig.  3 1  after  5PM  that  the  current  output  from  DER  increased  to 
almost  1  per  unit. 

D.  During  Resynchronization  to  the  Main  Grid  (at  8PM) 

At  8:00PM,  the  upstream  fault  was  cleared.  Fig.  32  illustrates  the  system  parameters  at  8PM. 

Islanded  mode  _  ^  Grid-connected  mode 
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Fig.  32.  Simulation  results  during  the  synchronization  of  the  microgrid  to  the  main  grid. 
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The  first  graph  in  Fig.  32  illustrates  the  fluctuation  in  voltage  and  current  at  the  microgrid 
interconnection  point  in  per  unit  from  0. 1  seconds  before  the  fault  was  cleared  and  0.2  seconds 
after  the  fault  was  cleared.  Fig.  32  also  displays  the  DER  output  power  in  kW  in  the  second 
graph;  and  in  the  third  graph,  the  voltage  difference  between  the  RMS  voltage  at  the  microgrid 
side  and  the  RMS  voltage  at  the  main  grid  side  in  volts. 

During  the  upstream  outage,  the  main  grid  voltage  was  collapsed,  while  the  microgrid 
maintained  its  voltage  at  1  per  unit.  Therefore,  the  voltage  difference  between  the  main  grid  and 
the  microgrid,  as  shown  in  Fig.  32,  was  240V  during  the  islanded  mode.  As  the  upstream  circuit 
was  restored  at  8PM,  the  main  grid  voltage  gradually  increased.  Therefore,  the  voltage 
difference,  as  shown  in  Fig.  32,  gradually  decreased  at  around  8PM. 

The  control  agent  monitored  this  voltage  difference,  as  well  as  the  frequency  of  both  systems. 
Once  it  detected  that  the  voltage  and  frequency  differences  remained  below  a  certain  threshold, 
the  control  agent  issued  three  control  signals.  These  control  signals  were  to  (i)  close  the  main 
switch;  (ii)  reconnect  non-critical  load  circuit  breakers;  and  (iii)  decrease  the  DER  output  to 
prevent  instantaneous  overshoot  in  the  system’s  voltage  and  current. 

In  this  case,  at  0.02  seconds  after  8PM,  the  control  agent  issued  the  synchronization  signal.  The 
system  voltage  was  again  stabilized  at  1  per  unit.  There  was  minor  current  overshooting  in  the 
first  half  a  cycle.  Then  again  the  system  was  stabilized.  The  DER  output  demonstrated  a  small 
ripple  for  about  0. 1  seconds  before  arriving  at  its  steady  state. 


C)  Case  study  3:  Secure  Selected  Non-Critical  Loads  during  Emergencies 

In  this  scenario,  it  was  assumed  that  the  outage  started  from  12:20AM  until  midnight  the  next 
day.  Since  there  was  excess  capacity  from  the  DER  unit  after  serving  critical  loads,  this  excess 
capacity  would  be  shared  by  some  of  the  non-critical  loads  in  the  community.  Here,  we 
presented  results  of  the  ID  APS  agent  management  system  that  allowed  the  limited  DER  capacity 
to  be  shared  among  non-critical  loads  during  outages. 

Before  the  outage  occurred  or  during  the  outage,  homeowners  or  facility  owners  were  asked  to 
enter  the  maximum  amount  of  money  that  ones  would  like  to  spend  on  non-critical  loads.  Users 
also  were  asked  to  select  the  time  of  the  day  that  they  would  like  to  operate  non-critical  loads. 
This  case  study  demonstrated  how  non-critical  loads  were  secured  within  the  ID  APS  microgrid, 
during  an  outage  that  lasted  almost  24  hours.  In  this  case,  we  focused  on  a  community  of  five 
houses.  Table  5  showed  the  $  limit  and  hours  for  non-critical  load  operation  selected  by  the 
owners  of  five  houses. 


Table  5.  Requirements  to  serve  non-critical  loads  during  outages 


House  1 

House  2 

House  3 

House  4 

House  5 

$  Limit 

$10 

$15 

$20 

$25 

$35 

Selected  non- 
critical  load 
operating  hours 

01:00-02:00 

05:00-10:00 

17:00-22:00 

01:00-02:00 

07:00-11:00 

19:00-22:00 

05:00-09:00 

17:00-21:00 

07:00-10:00 

19:00-22:00 

05:00-11:00 

17:00-23:00 
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As  a  standard  practice,  once  the  agents  detected  an  upstream  outage  at  12:20AM,  all  critical 
loads  were  secured.  The  decision  to  select  which  non-critical  loads  to  be  served  depended  on  the 
load  priority  (for  example,  non-critical  loads  in  critical  facility  like  hospital  are  of  higher 
priority),  amount  of  money  homeowners  were  willing  to  pay  for  their  non-critical  loads  to  be 
served,  and  their  desired  duration  of  operation.  Fig.  33  illustrates  additional  DER  output  to  serve 
non-critical  loads  when  needed. 
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Fig.  33.  DER  output  (kW)  to  serve  critical  and  non-critical  loads 
Fig.  34-38  illustrates  the  results  of  non-critical  load  operation  for  houses  no.  1-5,  respectively. 
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Fig.  34.  Non-critical  load 

■  Other 
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during  the  following 

■  Freezer 
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(a)  01:00-02:00 

(b)  05:00-10:00 

(c)  17:00-22:00 

Limitation  $10. 
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Fig.  35.  Non-critical  load 
operation  (house  2) 

In  this  case,  the 
homeowners  selected  to 
operate  non-critical  loads 
during  the  following 
hours: 

(a)  01:00-02:00 

(b)  07:00-11:00 

(c)  19:00-22:00 

Limitation  $15. 

Fig.  36.  Non-critical  load 
operation  (house  3) 

In  this  case,  the 
homeowners  selected  to 
operate  non-critical  loads 
during  the  following 
hours: 

(a)  05:00-09:00 

(b)  17:00-21:00 

Limitation  $20. 

Fig.  37.  Non-critical  load 
operation  (house  4) 

In  this  case,  the 
homeowners  selected  to 
operate  non-critical  loads 
during  the  following 
hours: 

(a)  07:00-10:00 

(b)  19:00-22:00 

Limitation  $25. 
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Fig.  38.  Non-critical  load 
operation  (house  5) 

In  this  case,  the 
homeowners  selected  to 
operate  non-critical  loads 
during  the  following 
hours: 

(a)  05:00-11:00 

(b)  17:00-23:00 

Limitation  $35. 


The  results  illustrate  that  the  ID  APS  energy  management  system  was  capable  of  controlling  DER 
and  loads  such  that  critical  loads  could  be  secured  during  emergency  and  some  of  the  non-critical 
loads  could  be  served  when  needed. 
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6.  Conclusions  and  Implications  for  Future  Research 


In  response  to  the  DoD’s  statement  of  need,  this  work  was  the  proof-of-concept  effort  to  study 
the  microgrid  operation  in  both  islanded  and  grid-connected  modes.  This  effort  resulted  in  the 
intelligent  distributed  autonomous  power  grid  that  could  facilitate  the  use  of  renewable  energy 
technologies,  thereby  minimizing  reliance  on  external  energy  resources  and  reducing  fossil  fuel 
consumption. 

The  following  five  tasks  as  proposed  were  accomplished: 

Task  1:  Development  of  ID  APS  device  models 

In  this  task,  models  of  ID  APS  physical  components  were  developed  in  Matlab/Simulink. 
These  models  included  solar  photovoltaic  modules,  wind  turbine  generators,  microturbines, 
fuel  cells  and  battery  storage. 

Task  2:  Development  of  local  control  algorithms  to  control  ID  APS  device  models 

Local  control  algorithms  were  developed  to  control  DER  units  for  their  grid  interface.  The 
results  indicated  that  the  control  unit  could  accommodate  the  integration  of  DER  units  into  the 
utility  grid.  Voltage  and  frequency  could  be  controlled  to  follow  those  of  the  grid  or  follow 
system  reference. 

Task  3:  Development  of  an  IDAPS  energy  management  system 

The  IDAPS  energy  management  system  was  developed  based  on  agent  technologies.  The 
IDAPS  energy  management  system  comprised  four  agents,  including  DER  agent,  user  agent, 
control  agent  and  database  agent.  Each  agent  had  its  objectives  and  responsibilities. 

Task  4:  Development  of  addressable  IP-based  communication  interfaces  for  connecting 
IDAPS  device  models  and  IDAPS  energy  management  system 

Under  task  4,  the  connection  between  the  IDAPS  physical  and  cyber  systems  was  successfully 
established.  The  simulation  conducted  under  task  5  integrated  this  feature  that  allowed  agents 
to  communicate  with  physical  elements  in  the  simulation  circuit  via  TCP/IP. 

Task  5:  Simulation  and  evaluation  of  the  microgrid  in  both  parallel  and  islanded  operations 

Several  scenarios  were  simulated  to  evaluate  the  IDAPS  microgrid  and  its  energy 
management  system  to  perform  demand  management  during  the  grid  connected  mode,  and 
secure  critical  loads  when  an  upstream  fault  was  detected. 


The  simulation  results  demonstrated  that  the  proposed  IDAPS  microgrid  allowed  various  DER 
sources,  including  renewables,  to  network  and  collaborate  in  an  efficient  manner.  This  work  met 
the  proposed  objective  by  illustrating  that  the  IDAPS  microgrid  could:  (i)  perform  demand 
management  during  nonnal  operating  condition;  (ii)  island  the  microgrid  from  the  main  grid 
once  an  upstream  fault  is  detected;  (iii)  secure  critical  loads  and  shed  non-critical  loads  according 
to  the  given  priority  list  during  emergencies;  and  (iv)  resynchronize  the  microgrid  to  the  main 
grid  after  an  upstream  fault  is  cleared.  Several  technical  challenges  were  addressed  by  this  work, 
including  analyzing  dynamic  response  of  agent-based  control  strategies,  defining  multi-agent 
system’s  role  and  definition,  and  evaluating  on-  and  off-grid  operation  of  a  microgrid  with 
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renewable  distributed  generation  systems.  Addressing  these  challenges  would  contribute  to 
facilitating  the  implementation  of  renewable  energy  sources  in  a  microgrid,  thus  enhancing 
energy  security  and  reliability  for  the  mission-critical  parts  of  military  bases  and  campus-type 
facilities. 

It  is  possible  to  extend  the  modeling  and  simulation  development  based  on  this  work  so  that  it 
can  be  used  to  explore  the  feasibility  of  the  development  of  microgrids  in  a  campus-type  facility. 
The  developed  simulation  platform  with  some  modifications  can  be  used  to  realistically  simulate 
and  validate  any  proposed  microgrid  generation  and  load  control  algorithms.  Therefore,  it  can  be 
used  to  quantify  the  risks  and  the  impact  of  implementing  the  proposed  control  algorithms  before 
the  actual  field  implementation.  Based  on  this  work,  several  microgrid  deployment  scenarios  can 
be  validated  and  guidelines  for  microgrid  development  can  be  provided. 
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Appendix  A  -  Publications 


This  work  has  resulted  in  one  Master’s  thesis  on  multi-agent  systems,  and  five  publications  in 
several  IEEE  conferences  with  full  paper  reviews. 

•  Feroze,  H.,  “Multi-Agent  Systems  in  Microgrids:  Design  &  Implementation”,  Master 
Thesis,  Virginia  Tech,  August  2009. 

•  Pipattanasomporn,  M.,  Feroze,  F.,  and  Rahman,  S.,  “Multi-agent  systems  in  a  distributed 
smart  grid:  design  and  implementation”,  In  Proc.  IEEE  Power  System  Conference  & 
Exposition,  March  18,  2009,  Seattle,  WA. 

•  Waseem,  I.,  Pipattanasomporn,  M.,  and  Rahman,  S.,  “Reliability  benefits  of  distributed 
generation  as  a  backup  source”,  In  Proc.  2009  IEEE  Power  &  Energy  Society  General 
Meeting,  July  2009,  Calgary,  Canada. 

•  Shao,  S.,  Pipattanasomporn,  M.,  and  Rahman,  S.,  “Challenges  of  PHEV  penetration  to 
the  residential  distribution  network”,  In  Proc.  2009  IEEE  Power  &  Energy  Society 
General  Meeting,  July  2009,  Calgary,  Canada. 

•  Zhang,  T.,  Shao,  S.,  Pipattanasomporn,  M.,  and  Rahman,  S.,  “Impact  of  variable 
electricity  rates  on  distribution  feeder  load  shapes  in  a  smart  grid”,  In  Proc.  IEEE 
PES/IAS  Conference  on  Sustainable  Alternative  Energy,  September  2009,  Spain. 

•  Shao,  S.,  Zhang,  T.,  Pipattanasomporn,  M.,  and  Rahman,  S.,  “Impact  of  TOU  rates  on 
distribution  load  shapes  in  a  smart  grid  with  PHEV  penetration”,  Accepted  for 
presentation  in  IEEE  PES  Transmission  &  Distribution  Conference  &  Exposition,  April 
2010,  New  Orleans,  Louisiana. 
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